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МАТЕМАТИКА 


УДК 517.951 


О НЕОБХОДИМЫХ И ДОСТАТОЧНЫХ УСЛОВИЯХ 
НЕТЕРОВОСТИ ДЛЯ ПОЛУЭЛЛИПТИЧЕСКОГО 
ОПЕРАТОРА СО СПЕЦИАЛЬНЫМИ 
КОЭФФИЦИЕНТАМИ 


А.А. Дарбинян, А.Г. Туманян 
Российско-Армянский университет 
аттаиЕг@уайоо.сот, атл4итапуаи @отай.сот 
АННОТАЦИЯ 


В данной статье исследована нетеровость полуэллиптического опе- 
ратора, с переменными коэффициентами, имеющими определенное поведе- 
ние на бесконечности. Доказана нетеровость полуэллиптических операторов 
со специальными коэффициентами в весовых соболевских пространствах. 

Ключевые слова: нетеровость, полуэллиптический оператор, ани- 


зотропные пространства. 


Данная работа посвящена исследованию нетеровости полу- 
эллиптических операторов с переменными коэффициентами в ани- 
зотропных пространствах Соболева. 

Исследованию нетеровости эллиптических операторов посвя- 


щена работа [1], где доказана эквивалентность нетеровости и эллип- 
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тичности в определенных пространствах на компактных многообра- 
зиях. В работе [2] в терминах предельных операторов получены 
условия для нетеровости эллиптических операторов в пространствах 
Соболева в ВП. 

В [3] доказана нетеровость полуэллиптического оператора с 
постоянными коэффициентами в анизотропных пространствах Собо- 
лева в ограниченной области. В [4] получено необходимое и до- 
статочное условие нетеровости для операторов с постоянными коэф- 
фициентами в анизотропных пространствах в”. Нетеровости 
полуэллиптического оператора в весовых пространствах посвящены 
работы [5—6]. 

Определение 1. Ограниченный линейный оператор А, опреде- 
ленный на всем банаховом пространстве Х и действующий в банахо- 
вом пространстве У, называется нетеровым, если выполняются сле- 
дующие условия: 

|) область значений оператора А замкнуто (1т(А) = Ип(А)); 

2) ядро оператора А является конечномерным (4йт. Кег(А) < о); 

3) коядро оператора А конечномерно (@йт сокет(А) = @туУ/ 
[т (А) < 09). 

Индексом нетерового оператора А называется разность между 
размерностью ядра и коядра: 

та(А) = атКег(А) — ат сокет(А). 

Рассмотрим дифференциальную форму: 

26:2) = ) в. (1) 
(1) $5 


гдез Е МиЕ М", (а: у) =“ +..-+ 0 Оо 


У1 


о: 10 
[ о = (х.,..,х,) Е В". Обозначим через 
Е 


В(ь0)= ) аабдр“ (2) 


(ам)=5 
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главную часть дифференциальной формы Р(х, 2), а через 
оз(, 8) = У азбдЕа (3) 
(а/)=5 


символ Р.(х, О). 


Определение 2. Говорят, что дифференциальная форма 
Р(х, 2) полуэллиптична в точке ху Е В", если 0.(ху,&) = 0 при всех 
ЕВ", || = 0. 

Определение 3. Говорят, что дифференциальная форма 
Р(х, 0) полуэллиптична в К", если Р(х,0) полуэллиптична в каждой 
точке х Е АП. 

Для Ае №,уеЕ М” обозначим: 


СКУ := СЕ(Вп) 
и [@6о: РВа(х) Е С(В"), зир |[РЁа(| < ©, УВ 
ХЕВП 


Е п, (В: < к 


О” — 0” (в”) 
_ ‚ |296) 


о 
ра Ки. п 0. 
| [бо ес :9(00) > 0,УхЕ В ао 0 о 


3 Опри |х| > ® УВ Е 7#, (В: м) <«| 


Дляу Е М, КЕМиаЕ 0*"через Н”” и НЕ”обозначим, соот- 
ветственно, множества измеримых функций {и} с конечными нор- 


мами: 


Пи, = У ИРби све, 
(а:у)<К 


[р = о |2“ Г о 
(ам) <К 
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Условие Т. Пусть КЕ М, К > $ и коэффициенты аз(х) главной 
части дифференциальной формы Р(х,0) при (а:у) = $ постоянны, 
ав(х) Е С*` 8” при (а: у) < $ и, при этом, существуют постоянные йд 
такие, что РВ(а„(х) — аи ) 3 0 при |х| -> © для всех В Е 71, (В:у) < 
К 5. 

Пусть дифференциальная форма Р(х,0) вида (1) удовлетворяет 
условию Т Тогда Р(х,О) порождает ограниченный линейный 
оператор, действующий из Н*” в НА-5" (Р:Н®" > НК 5”). 


Для дифференциальной формы Р(х,0), удовлетворяющей ус- 


Р(р) = р ав0® + о а рй. 
(ам)=5 (а:у)<5 


Для МЕВ.,С ЕЁ" и УЕМ" обозначим Км = {хЕ В": [х| < 


1 
УЕ = 0.8"). 

Теорема 1. Пусть дифференциальная форма Р(х,0) удовлет- 
воряет условию 1. Тогда Р(х,0): Н®" > НЕ” является нетеровым 
тогда и только тогда, когда существует положительная постоян- 
ная д такая, что 

|Р(2)| > 61+ 18°, УЕ Е В". 

Доказательство. Обозначим 

т,0) = > (а. бд —а,)0“. 
(а:1/)<5 
Имеет место следующее представление: 
Р(х,0) =Р(Рр) + Т(х,О). 

Докажем компактность оператора Т(х, 2): Н®" — НК», 

Так как ОВ(а(х) — аи) 8 0 при |х| — 0 для всех (В:у) < К - 5, 
(а: у) < 5, то для произвольного = > 0 существует №(=) > 0 такое, что 

[28 (а›(х) — а. )| < &,Мх Е В"\Кицьу, (В: У) < К - 5, (в: у) < $. 
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Пусть ф(х) Е Су’ такая, что 0< ф(х) < 1иФ(х) = 1 при [х| < 1 
и Ф(х) =0 при |х|>2. Тогда ф:(х):=Фф (5) Е Су’, 5иррф: с 
Ком(=)» Ф=(х) = 1 прих Е Кк(;). 


Обозначим Т; = (1-Ф2Т, Т/ = Ф.Т. Тогда с некоторой по- 

стоянной С > 0 

Тем кз» = И — ФТИ < СЕЙ, Мм Е НАУ, 
аТ.’ = ф-Т - компактный оператор из Н*” в Н*`5”. Следовательно, в 
силу теоремы 8.3.2 работы [7] оператор Т(х, В): Н®" > Н*7” ком- 
пактен. 

Поэтому на основании теоремы 8.5.10 работы [7], оператор 
Р(х, 0): Н®” -› Н®-5” является нетеровым тогда и только тогда, когда 
Р(Р): Н®" > Н*-5” являетя нетеровым. 

В меореме 3 работы [4] установлено, что оператор Р(р): Н*” 
НЕ нетеров тогда и только тогда, когда с некоторой постоянной 
5>0 

|Р(5)|> ба +8, УЕ Е В". 

Замечание 1. Так как из условия [2(2)| >0(1+12|,)°, У Е В" 
следует полуэллиптичность Р(2), то из нетеровости Р(х, 2): Н®" 
НЮ следует, что Р(х, О) также полуэллиптичен в ВП. 

Условие П. Пусть К Е М,К > 5,4 Е О*” и коэффициенты аа(х) 
главной части дифференциальной формы Р(х,0) при (а:у) = 
$ постоянны и при (а: и) < $ имеет место следующее представление: 

аа (о) = а0(2)960°-@) + а! о, 
где а0(х) и аЁ (х) удовлетворяют следующим условиям: 
1) 28 (а1 (х)) = о(а (к) @м +") при |[х| — 9 для всех 
(а: у) <5и(В:у) <К-$; 
2) а°(х) Е С" 5" и существуют @х, константы такие, что а0 (х) 3 


а при |х| — со для всех (а: т) < $. 
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При условии П дифференциальная форма Р(х,))) порождает 
ограниченный линейный оператор из На "вН о зд 
Для дифференциальной формы Р(х, 0), коэффициенты которой 
удовлетворяют условию П, обозначим 
Р(х,0) = а а„0“ + о п РОЯ: 
(а:у)=5 (а) <$ 
Условие Ш. Пусть 
о а_а + о Я.А) а О, МЕ Е В", > 0. 
(а) =5 (а:)<$ 
Теорема 2. Пусть дифференциальная форма Р(х,2), коэффи- 
циенты которой удовлетворяют условиям П, Ш, полуэллиптична в А". 
Тогда Р(х, О): я." > Я ” нетеров. 
Доказательство. На основании условия П представим Р(х,0)в 
следующем виде: 
Р(х,0) = Рбк,0)+ > (ад - даб ря + У  абдр“ 
(ал) <$ (а:/)<5 


и обозначим 


Т1(х,0) = й (а (©) — а )а60*-@^ ра, 


(а) <$5 
Т?(х,0) = р ат (х) р“. 
(а:/)<$ 
Так как для всех (@: у) < 5а°(х) Зави 30 при |х| — 00, то 


90) 
для произвольного = > 0 существует № (Е) > 0 такое, что 
127 (а<(х) — @в)| 
909") 
Пусть Ф(х) Е С5° такая, что 0 < ф(х) < 1иф(х) = 1 при [х| < 1 


< & УХЕ А"\Км (=), (у: м) < К — $, (у) < 5. (4) 


х 
№1 (Е) 


Км, (=) ф(х) = 1 прихЕ Км, (=) и! = те ФТ 55 и. Г. 





и ф(х) =0 при |х| >2. Тогда ф.(®):= Ф( ) Е С0°, зирро; с 
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В силу оценки (4) с некоторыми постоянными С.,С› > 0 для 
всех (р:у) < К—$и (а: м) < $ имеем 
[28 ((1-— Ф-©0) (або - а) “ис а | 


128") 
<) |2” ((1- 9-62) 
0=у<В 
о) оао | 
2 








мт. Е ина о: ее 
зо т [2°+8-Ум | 
0<у<В 
Отсюда с некоторыми постоянными (5, С. > 0 получим, что 


[|1 =: ф=) Те х Ви к-за 
< >) У |8 ((@ - бд) (ад 


(Ву) =К-$ (а1/)<$ 


— би) р“ и(х)) 400"—*_@ 


12(В")" 








< С.=м 
ео 5 СЕ 


\Ми;ЕН*. 

Так как Ри НА "> Но и. компактен, то отсюда на основании 
теоремы 8.3.2 работы [7] получим, что оператор Т*(х,Б): Н* > 
Н "также компактен. 

Рассмотрим оператор Т? (х, 2): В > И — 

В силу пункта 1) условия П для произвольного = > 0 существу- 
ет №. (=) > 0 такое, что 

12"а&(х)| 
або ато) < =, Ух Е В"\Км, (5), (У:У) < Е — $, (а: ) < 5. 


Тогда 
Феб = ( 
Км, (=). 


х 
№2 (=) 





) Е Со’, 5ирру: с К>м. (=), ф:(х) =1при хЕ 


Проводя аналогичные рассуждения как при доказательстве ком- 


пактности оператора Т\(х, В): Но" > И с заменой Т\ на Т? и Ф; 
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на ., в силу теоремы 8.3.2 работы [7] получим, что Т?(х, В): в. > 


К-5У 
На 


так же компактен. 


Так как при уловиях теоремы в силу теоремы 2 работы [5] опе- 


ратор Р(х, В): Ни" > НЕ” нетеров, то на основании теоремы 8.5.10 


работы [7] получаем, что оператор Р(х, 0): На” эн: ® так же 


нетеров. 


Работа выполнена в рамках тематического финансирования РАУ 


из средств МОБН РФ и при поддержке тематического финансиро- 


вания комитета науки при МОНРА (код проекта 5С$ №15Т-1А197). 
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АЛГОРИТМ СЛЕЖЕНИЯ ЗА ТРАССОЙ 
ПО ВИДЕОСЪЕМКАМ БЕСПИЛОТНОГО 
ЛЕТАЮЩЕГО АППАРАТА 


Д.Г. Асатрян, С.М. Овсепян, Т.Б. Айрапетян 


Российско-Армянский университет 
АННОТАЦИЯ 


В настоящее время беспилотные летающие аппараты (БЛА) находят 
применение во многих сферах деятельности человека. Обнаружение и сле- 
жение за трассой при помощи видеосъемки БЛА играют важную роль в раз- 
личных задачах мониторинга местности, автоматического обнаружения 
объектов интереса и др. В настоящей работе исследована работоспособность 
предложенной ранее процедуры слежения за трассой, используя передовые 
методы обработки изображений, основанные на использовании структурных 
свойств изображений. Процедура основана на многократном последователь- 
ном применении алгоритмов двух типов — алгоритмов оценивания сходства 
любых участков изображения и их доминантного направления. Работоспособ- 
ность предложенной процедуры иллюстрирована численными примерами. 

Ключевые слова: Беспилотный летающий аппарат, видеосъемка, гра- 


диентное поле, сходство изображений, доминантное направление изображения. 


1. Введение 

В настоящее время наблюдается бурный рост приложений бес- 
пилотных летающих аппаратов (БЛА) в самых разных сферах деятель- 
ности человека. Применение в составе БЛА современной навигацион- 
ной, видео- и коммуникационной аппаратуры делает БЛА незамени- 


мым инструментом для решения многих задач. 
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Одним из важных применений БЛА является их использование 
для мониторинга дороги, контроля траффика, слежения находящихся 
на ней объектов, автоматической навигации и др. Слежение за трассой 
является автоматическим или полуавтоматическим процессом, кото- 
рый использует видеоинформацию, получаемую БЛА и специальные 
методы обработки изображений. 

Сложность этой задачи связана с наличием многочисленных 
факторов, влияющих на качество и содержание изображений на кад- 
рах видеосъемки. Основными факторами, постоянно действующими 
на процесс снятия видеоизображения, являются эволюции аппарата, 
турбулентность атмосферы, наличие ветров, вибрации борта аппарата, 
а также внешние преднамеренные и непреднамеренные помехи и иска- 
жения изображения реальной сцены. Все это приводит к тому, что в 
соседних кадрах видеосъемки изображение одной и той же сцены, рав- 
но как и его геометрические параметры, отличаются непредсказуемым 
образом. 

Для эффективной обработки последовательности изображений, 
искаженных вследствие воздействия описанных выше факторов, необ- 
ходимо разработать специальные алгоритмы, достаточно устойчивые 
относительно подобных искажений. Как известно, такими алгоритма- 
ми пользуется человек, рассматривающий последовательность кадров 
видеосъемки, поэтому при разработке эффективных алгоритмов авто- 
матической обработки желательно использовать соответствующие 
свойства визульной системы человека (ВСЧ). 

Следует заметить, что использование глобальных систем пози- 
ционирования СР$ или СГОМА$$ и других типов навигационной тех- 
ники, установленной на БЛА, существенно облегчает решение расс- 
матриваемой задачи, однако в последние годы становятся актуальны- 
ми задачи, в которых использование этих средств нецелесообразно 


(для обеспечения скрытности) или невозможно (при выходе из строя 
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по разным причинам соответствующего оборудования). В подобных 
случаях важной и чуть ли не единственной информацией для выполне- 
ния решаемой задачи является видеосъемка местности и он-лайн обра- 
ботка последовательности видеоизображений. 

В настоящей работе сделана попытка разработки методики сле- 
жения за дорогой на основании имеющейся при полете БЛА видеоин- 
формации, дополнительно используя лишь информацию электронного 
компаса для определения направления движения. 

В литературе предложено много подходов и методов, относя- 
щихся к задаче поиска, обнаружения и слежения за дорогой и/или на- 
ходящихся на ней объектов [1-3]. При этом используются различные 
признаки, характеризующие присутствие на изображении участков до- 
роги или объектов на них, такие как движение [4], фон [5—6] и др. 

Мы ограничиваемся рассмотрением ситуации, когда начальный 
участок дороги уже найден или задан, поэтому в настоящей работе 
вопрос поиска и обнаружения дороги в общем случае не рассматри- 
вается. Вместо этого предполагается, что в начале поиска имеется об- 
разец, изображающий какую-нибудь характерную часть дороги и ее 
окрестность, а в процессе поиска в каждый момент времени видеока- 
мере БЛА доступен только небольшой участок дороги, содержащийся 
в текущем кадре. Следовательно, процедура слежения за дорогой сво- 
дится лишь к использованию информации о свойствах дороги, содер- 


жащейся в упомянутом выше образце. 


2. Описание методики обработки и результаты 
экспериментов 
Как отмечено выше, методика не предусматривает использование 
глобальных систем позиционирования ОРЗ или ОТГОМА$$, следова- 


тельно, решение принимается по результатам обработки видеопоследо- 
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вательности, а корректировка курса полета производится с учетом пока- 
заний электронного компаса, установленного на борту БЛА. 

Процедура обработки последовательности кадров, доставленных 
камерой БЛА, основана на многократном применении двух алгорит- 
мов — алгоритма оценивания сходства и алгоритма оценивания доми- 
нантного направления изображений рассматриваемых сцен. Эти алго- 
ритмы предложены ранее в [10-12] и показали эффективность их при- 
менения в различных приложениях. Оба алгоритма основаны на мате- 
матических моделях, использующих различные характеристики гра- 
диентного поля изображения. В [13] описана схема использования 
этих алгоритмов в задаче слежения за дорогой. Однако обоснование 
работоспособности предложенной в этой работе процедуры на основа- 
нии численных расчетов не приведено. 

Вначале кратко напомним суть упомянутых двух основных алго- 
ритмов. 

Алгоритм оценивания сходства двух изображений использует 
матрицы магнитуд градиентов, полученных применением оператора 
Собела к интенсивностям пикселов, принимая элементы матриц за 
случайную величину, распределенную двухпараметрическим законом 
Вейбулла. Мера сходства при этом рассчитывается по статистическим 


оценкам параметров (п, с ) этих распределений. Соответствующая ме- 
2 =. 
ра сходства \/" определяется формулой: 


р’з_ ПИКАП, тай.) 


= 0<И” <1. 
тах(И, ‚|, ) тах(о,,с,) 


Алгоритм оценивания доминантного направления изображений 
основан на построении эллипса рассеяния градиентного поля, опреде- 
ляемого по матрицам вертикального и горизонтального градиентов 
изображения. Направление доминантного направления изображения 
отождествляется с направлением главной оси эллипса, задаваемого 


уравнением: 
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2 
= (и [Ин _ Зри (вн Ниву Ши), (вм) в 
1 Рау Он ОнОбу 5. 
где Цин г | , Он р Оу — выборочные средние и среднеквадратические 


отклонения горизонтальной И вертикальной составляющих градиент- 


ного поля изображения, р‚„„ — коэффициент корреляции между ними, 


С- константа. 

В [13] приведены как расчетные формулы, так и некоторые чис- 
ленные примеры указанных двух характеристик. 

В настоящей работе на конкретном примере приводятся резуль- 
таты обработки различных участков изображения, содержащих пре- 
имущественно дорогу или фон (в данном случае — лесной массив), по- 


казывающие возможность БЛА выполнять необходимые маневры для 


слежения за дорогой на основании этой информации. 





Рис. 1. Изображение дороги. Прямоугольники обозначают участки, 
содержащие преимущественно дорогу (пронумерованы цифрами 
красного цвета) или фон (пронумерованы цифрами белого цвета). 


На Рис. 1 приведено изображение местности с дорогой, а также 


выбранные на нем участки, содержащие преимущественно часть доро- 
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ги (пронумерованы цифрами красного цвета) или фона (пронумерова- 
ны цифрами белого цвета). Размеры участков примерно одинаковы, но 
не фиксированы (имитация реальной ситуации при полете БЛА). 
Смысл этой иллюстрации в том, что БЛА в последовательные момен- 
ты времени реального полета может находиться над участком, имею- 
щем вид типа одного из изображенных в прямоугольниках. В каждый 
момент полета решение о выборе направления дальнейшего движения 
БЛА может принимать в зависимости от значений рассматриваемых 
двух характеристик — степени сходства текущего участка с заранее 
выбранным участком (обозначено номером 1) и направлением доми- 
нантного направления изображения текущего участка (угол между до- 
минантным направлением относительно горизонтальной оси). Эти зна- 
чения приведены в Таблице. В левой половине Таблицы размещены 
данные, относящиеся к участкам с дорогой (участки с номерами от 1 до 
7), а в правой размещены участки с фоном (с номерами от 8 до 14). 

Таким образом, предложенная методика представляется вполне 
работоспособным и может быть реализована в системах программного 
обеспечения БЛА при решении задачи слежения за дорогой с исполь- 
зованием лишь результатов видеосъемки и показаний электронного 
компаса без привлечения навигационных систем. 

Сравнение данных Таблицы показывает, что значения \” из ле- 
вой части Таблицы существенно больше значений из правой части. 
Это различие позволяет при корректировке направления движения 
БЛА с уверенностью выбрать решение в сторону участка с большим 
значением \№ и, таким образом, осуществлять слежение за дорогой 


(см. также алгоритм корректировки, описанный в [13]). 
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Таблица. 


Значения меры сходства выбранных участков с участком №1 


и доминантное направление участков (угол в градусах). 














№ \? Угол, град. | № участка \? Угол, град. 
участка 
1 1.00 5.7 8 0.45 163.0 
2 0.91 11.56 9 0.42 164.1 
3 0.89 16.75 10 0.59 177.2 
4 0.81 24.03 И 0.46 174.9 
5 0.81 29.33 12 0.51 152.9 
6 0.91 35.95 13 0.47 164.5 
7 0.86 42.42 14 0.51 169.8 























Аналогичные выводы можно сделать, анализируя значения до- 
минантного направления участков: первые 7 участков с большими зна- 
чениями \\"? показывают плавное изменение доминирующего направ- 
ления, что вполне согласуется с соответствующим результатом ви- 
зуального анализа. На Рис. 2 показано изображение с выбранными 
участками, на которые нанесены стрелки красного цвета с углом нак- 
лона, взятом из Таблицы. Поскольку при выборе направления движе- 
ния БЛА углы доминантных направлений с низкими значеними Ув 
расчет браться не будут, на этом рисунке стрелками белого цвета по- 


казаны направления движения БЛА после корректировки. 
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Рис. 2. Иллюстрация расчетного направления движения 


по информации Таблицы. 


3. Заключение 

Решается задача слежения за дорогой, используя лишь инфор- 
мацию по видеосъемке БЛА местности с дорогой и показания элект- 
ронного компаса. Показана возможность решения этой задачи путем 
последовательного многократного применения процедур оценивания 
сходства и доминантного направления различных участков изображе- 
ния. Приведенные численные расчеты показывают работоспособность 


предложенного подхода. 
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оп фе зиисфга! ргорегиез оЁ ппасе. ТВе ргоседиге 15 Базе оп шире апа 
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24 Алгоритм слежения за трассой по видеосъемкам беспилотного летающего аппарата 





ПоОПОЭЛЬ ©0201, ОЧГРЬ $6 ОПУЧИОГОХОЗОО ОБЭПЗПЧ, 
СОЗГПБЛЛЬ ПРАВУ ОТаПГБоО 


2..9. Чицлорлае, 9.С. 2оЧаб фу, $.Р. Дацрана вилл 
ЦОФПФПЬО 


Збрбицот а шмота рврот ларрёро (099) ЧррлаоЧотч Ба Чарз- 
Чинре дорботабоелай рааЧиаер поро броы: 0.90 шбамаарлаьо- 
Чл одботераИр Чароттот Боллобар6 ра И ооо ЧЕроцебрр Чилр- 
бор обр Ба ато шблымарр Баба, БЕиаррорронешмт а Прада 
ОРЮР оо Чао Боллобарер ола 6 о бора роту: Осца шориа- 
иа@ро Ббиаари ЧЕ Е Шарочот ола ` бора Ораоьи чаролоар аб 
ара арар чарами ооо р, одимааорбшю] ШаротоЧкабраира 
Боло ое юри офилааорбЫшт Чрул Бр чб` имо Еробрр Чаи 
пара ор ЧбЕорщабр: Сормаоа ара БУ Чиаб Е рот шар ччаорррч- 
ЧЫрЬ Чрубоа о | Баопрца ла ИБроло ле Чрал, оробр 6&` манор рр Ца 
Чл ло брр ЧИаботерай © имо рр 9629200 пыпоторлаа да- 
Боло ааа ицарорре( бр: ДомомарЧилб рава дарар 12 ааморот алое - 
пр оптом Е рЧачри опр роч: 

2Врио бр оботдиутт. розпт ширр, аи лари льу, ара - 
пара голом, прало О ЕрОЕрЬ Убаонарота, моло рр 96092001 пыдот- 
от: 


Вестник РАУ № 2, 20176, 25-40 25 





УДК 517.968.4 


О РАЗРЕШИМОСТИ ОДНОГО КЛАССА 
НЕЛИНЕЙНЫХ ИНТЕГРАЛЬНЫХ УРАВНЕНИЙ 


НА ВСЕЙ ПРЯМОЙ” 


Х.А. Хачатрян', А.А. Сисакян” 


1 
Институт математики НАН Армении 
7 й ь ь 
Армянский Национальный Аграрный Университет 


Красй82@хтатЫег.ги, 55акКуапб@тай.ги 


АННОТАЦИЯ 


В настоящей статье исследован один класс интегральных уравнений на 
всей прямой со степенной нелинейностью. Указанный класс уравнений 
имеет применение в кинетической теории газов. Доказано существование 
нетривиальных и ограниченных решений. В конце работы приведены част- 
ные примеры указанных ядер. 

Библиография: 12 наименований. 

Ключевые слова: факторизация, последовательные приближения, ог- 


раниченное решение, монотонность, ядро. 


$ 1. Введение 
Настоящая работа посвящена изучению следующего класса нели- 


нейных интегральных уравнений на всей прямой: 


1*(х)= | ко) (в) хеВ (1) 





° Работа выполнена при финансовой поддержке ГКН МОН РА в рамках научного 
проекта № 5С$ 167Е-1А002. 
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относительно искомой непрерывной и вещественной функции / (х) : 


В уравнении (1) р>2 является нечетным числом, а ядерная функция 


допускает следующее представление: 
Ь 
К (х,1) = ‚а (хе “Чо ($), (х,)ЕВхХВ, (2) 


где о ($) — монотонно неубывающая на [а,5) функция, а функция 


&(х,5) определена на множестве Вх [а,Б) (0 <а<ь< оо) и 
удовлетворяет следующим условиям: 

А) аеС(Вх(а,Ь)), 

В) существует такое число 0 >> 0, что 


а(х,5) > 6, \(х,5) Е Вх[а,Б), зар а(х,5)<20, (3) 


(^у)ЕВХа,5) 
С) а“ (х,5 ) симметрична по первому аргументу: 
а(—х,5)=а(х,5), х> 0, 5Е |,ь), (4) 
О) существует число Г>0, являющееся основным периодом для 
функции &(х,5) по аргументу Х ‚ т.е. 
а(х+Т,5) =а(х,5), У (х,5) Е Вх (а,Б). (5) 
Решение уравнения (1) мы будем искать в следующем классе 


непрерывных на № функций: 
И= [Ф ф(-х) = —Ф(х), х>0, фЕ С(Е) ПД, (®)}}, (6) 
где РЁ» (В) — это пространство существенно ограниченных на Ю 


функций. 
Уравнение (1) возникает в кинетической теории газов (см. [1]). В 


линейном случае, когда р=|1 и К (х,1)=0 при #<0, вопросы 


разрешимости уравнения (1) обсуждались в работах [3, 4]. 
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В случае, когда Я (х,5) = СОП51, легко можно убедиться, что 
уравнение (1) сводится к линейному дифференциальному уравнению 
Риккати. Отметим также, что в случае, когда ядро К уравнения (1) 
имеет вид: 


К (х,г) ЕЕ Е. : (7) 


Ул 
то уравнение (1) с граничными условиями 
У (+) =+1 (8) 
описывает поле тахионов для открытой р-адической струны (см. [5, 
6]). Вопросам разрешимости граничной задачи (1), (7), (8) посвящены 
работы академика В.С. Владимирова (см. [6-9] и ссылки в них). Нако- 


нец, когда ядро К (х,г) зависит от разности аргументов, имеет об- 


щую структуру и удовлетворяет условиям 
+ со 
К (т)> 0, тЕК, /[ к(т)ат=ь 


+ оо 
К ЕЁ. (В), [т°к(т) < +со, К (—х) = К(х), ЕВ [0, оо), 
то одним из авторов настоящей работы доказана теорема существова- 
ния и единственности граничной задачи (1), (8) (см. [10]). 
В данной работе при условиях (2), А)-С), доказывается существо- 


вание нетривиального ограниченного и нечетного решения уравнения 
(1). В конце работы приведены частные примеры функций @ (х,5) , 
удовлетворяющих условиям А)-О). При выполнении дополнительного 


условия 0) устанавливается существование однопараметрического 


семейства нетривиальных решений. 
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$ 2. Обозначения и вспомогательные факты 


2.1. Сведение уравнения (1) к интегральному уравнению 

на полуоси 

Сперва уравнение (1) приведем к уравнению на положительной 
полупрямой следующим образом: пусть А (х) — непрерывное реше- 


ние следующего нелинейного интегрального уравнения на полуоси: 


Ех) = [(Кьй-КИы) Е а, хев", = © 


где К 68) задается согласно формуле (2), а Ку, (х,г) имеет следую- 


ЩИЙ ВИД: 
(6) ›= х,з)е “+ 405), (хбЕВхВ-. 

Тогда а 
Х(х)= 


является непрерывным решением уравнения (1), где Ё решение 


Е(х), если х> 0, 
(10) 
Е(-х), если х < 0, 


уравнения (9). 


Действительно, в силу условий А) и С) из уравнения (1) будем 


при 0 
[код а= | ко) г(а+ | койг(оа= 


= [к(х,-0/(-да+ | кодЕ( а | (К(хи)- Ки) Е (да = 


0 0 0 
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если же х<0 (-х>0), то в силу условия С) и определения Г 


имеем: 


ый Е ба [к —0/(-а м >) Е(@- 
кодека = [(кзы)- к) (9 = р) 


В $3 (пункт 3.1) мы докажем, что Ё(0)=0. 


2.2. Об одном вспомогательном уравнении Вольтерра 
Рассмотрим следующее линейное неоднородное интегральное 
уравнение Вольтерра: 


со 


4х) = &(х)+ [Уса хЕВ* (11) 


относительно искомой функции 4 (х) ‚ где 
Ь 
29 | г“ 4о(5), хеВ*, (12) 


У(х,г)=2К(х,1)0(1—х), У(хИЕВ ХВ", = (13) 
а 9- известная функция Хевисайда. 
В дальнейшем будем предполагать, что монотонная функция О’ 


удовлетворяет следующему условию: 


Ь 
а =—. 14 
[49(5)=5 (14) 


Уравнение (11) перепишем в операторной форме: 


30 О разрешимости одного класса нелинейных интегральных уравнений на всей прямой 





(1-У)ф= &, (15) 


где [ — единичный оператор, а Г— интегральный оператор Вольтерра 


следующего вида: 
= [Убол(0аь хЕК*. (16) 


Можно проверить, что оператор Г действует в каждом из про- 


странств В. где Е — одно из следующих банаховых пространств: 
Е, (В+), 15 р<-+оо, 1.(В*), С(В*), С,(В*), 
где С. (®-) — пространство непрерывных и ограниченных функций 


на В’, а © (®-) — пространство непрерывных функций на ®\, 
имеющих нулевой предел в +. 

Это следует из свойств А)-С), накладываемых на функцию 
а(х,5). 

Рассмотрим следующую задачу факторизации: для оператора И 


найти такой интегральный оператор И’ вида: 


5) Ре (х,1) (Г) 4 хЕВ', ГЕЕ, (17) 
Ш (х,1 ) ЧЕ < + оо, (18) 
х>0 
чтобы имело место Е 
1-Г=(1-#)(1-0), (19) 
где 
98/2 9 (Па, хеВ+, ГЕЕ. (0) 


Согласно результатам работы [11] имеем: 
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У? (х, 1) )2 х,5) Ве “40 (5 )0(#—х), 
У (х,1)Е "ЖК". (21) 
Из представления (21) и условия (3) о что 
— 
= У? (х,г) < | а 22 
нм убыаа [а)-з. |. 
р, ЛХ (х,#) раз(а х,5)- В) 40 (5 д < =1. (23) 
10 В 
Используя факторизацию (19), уравнение (15) запишем в следую- 
щем виде: 
(7-И)И-И)4= 5. 4) 
Обозначим 
(1-()9=4. (25) 
Тогда 
ЕЙ ав (26) 


или в раскрытом виде: 


(= (4) + Гуьда (да, хев*. о 


Так как ру < р <Ти ЕЁ (в) пд. (®-) ‚ то уравнение 








(27) имеет единственное решение в /, (в) ПД. (в ). 


Теперь займемся решением уравнения (25). Решение уравнения 
(25) задается посредством следующей формулы: 


4(х)=а(х) + 8 [а(#)4ь хеВ*. (28) 
Покажем, что 


т, (4) = ‚[ ха (хак < +. (29) 
0 
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Сначала убедимся, что 7, ( =) < о. Действительно так как 9>0, 


то из определения © (см. (12)), имеем 


(5)= [х( ке 2х: (5) < +. 
0 а0 
С другой стороны 


Следовательно, 


и х,1) 2 дама [в (9) 1) < +. 
00 


(30) 


Тогда, так как т, ( &) < +=, то из (27) в силу (30) имеем 


т, (4) < +=. 


Используя это, получим, что уравнение (25) имеет решение из 


пространства Г, (®*) Ра (в* ) ‚ме 
иЕЁ(В*) пд. (В*). 


2.3. Построение нетривиального решения однородного 
уравнения Вольтерра 


Рассмотрим следующее однородное уравнение Вольтерра: 


5%) = [6-5 а, хе, 


где 
Ь 


(жк) =2 [а(жз)е “о ($, (хи) ЕВ ХВ". 


Так как при (х1)ЕВ*хВ*, < х+Е, то 


(30 


(32) 


(33) 
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К(х1)> К,(х,ё), У(хр) ЕВ ХВ*. (34) 
Итак, неравенство (34) установлено. Из (34) следует, что 
у(х,1) > $(х,1), У(х) ЕВ ХЕ. (35) 


Следующая теорема в наших дальнейших рассуждениях сущест- 
венно понадобится. 
Теорема 1. При условиях А)-С) и (14) уравнение (32) обладает 


неотрицательным нетривиальным и существенно ограниченным ре- 
+ 
шением 5, причем 1-5 ЕЁ (В ) 


Доказательство: Рассмотрим следующее неоднородное уравне- 
ние Вольтерра: 


0(%)= («+ [бд -ъдоа, хев. = в9 


Для уравнения (36) введем следующие последовательные приближе- 
НИЯ: 


О, (х)= 8(х)+ | (5) -(х1)0, (4, хе", 67) 


(ХЕ н= од 


Индукцией по И нетрудно убедиться, что 


























О, (х)Т поп, хЕ к, (38) 
0, (х)<и(®) , п=0,1,2,.-., хЕВ*, (39) 
функции О, (х) измеримы по Хх. (40) 


Следовательно, последовательность функций ТО, (х)} о Имеет 
и= 


поточечный предел, когда И — со: Пт О, (х) =О(х), причем 
й—>=> 


#(х)< О(х)<и(х), хЕВ*. (41) 
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Согласно предельной теоремы Б. Леви (см.[12]), функция О\(х) 
является решением уравнения (36). Так как ИЕ [/ (^*) ВР (^* ) ‚ то 
из (41) следует, что ОЕ [ (®*) АЕ (^* ). 


Прямой проверкой можно убедиться, что уравнению (36) 














удовлетворяет также (), (х) =]. Индукцией по И легко можно 
доказать, что 
О, (х) < п=0,1,2,.-., хЕВ®. (42) 
Следовательно, 
0(х)< 0.(х), хЕВ* (43) 


Заметим, что 

0(х) О. (х), (44) 
ибо ОЕД(В*) ПЕ. (В*), а О. (х 
5(х)=0,(х)-О(х)=1-0(х)> 


является решением уравнения (32), причем 0<5 (х) <1: 5 (х) Е 0. 


) =1. Очевидно, что 


Теорема доказана. 
$3. Разрешимость уравнения (1). Примеры 


3.1 Построение нетривиального и ограниченного решения 
уравнения (9) 


Рассмотрим следующие последовательные приближения для 
уравнения (9): 


Её (х)= | (К(х)— Ко) Е(аь хе 45) 
0 


0 (х) =1, и=1,2..... 
Индукцией по И сначала убедимся, что 


Е, (х){ поп, хЕВ*. (46) 
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Сперва докажем, что РУ (х) <Е\(х), хЕВ". С учетом (2), (14) и (34) 


из (45) имеем 


Е” (х) = (Код к, бы) < [Кора < 


отсюда следует, что В (х) <1=А (х). Предполагая, что 


Е (х)< Е, (х) при некотором натуральном И и учитывая не- 


четность числа р >2, из (45) получим 


Е (х)< (Ко) - Кб) Е = ЕР (о), 


0 
откуда следует, что 


Теперь докажем, что 
1 


Та 
в6)>|5] 5(х), п = 0,1,2...., хЕВ”, (47) 
где 5(х) — решение однородного уравнения Вольтерра (32). 


Выполнение неравенства (47) в случае И = 0 сразу следует из (45) и из 
двойного неравенства 


0<5(х) < хЕВ" (48) 


(см. доказательство теоремы 1). Предположим, что (47) выполняется 


при некотором ИЕ М. Тогда в силу (34), (32) и (48) из (45) имеем: 
1 


в [5] (ко) -Кде)5да= 


х 


Р 


Поскольку р>2 - нечетное число, то из полученного неравенства 
следует: 
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1 
1 р 
(х)>|-| 5(х). 
Итак, (47) доказано. 
Используя непрерывность функции &(х,5) (см. условие А)), 


индукцией по И легко можно убедиться, что 
РеС(В*), и (49) 














Таким образом, из (46), (47) и (49) заключаем, что последовательность 
непрерывных функций [Е, (х)! ‚ Имеет поточечный предел, когда 
и= 


п—> =: ШиР, (х)=Е(х), причем 


5] 5(х)< Е(х)<1, хЕВ*. (50) 


Так кк РЕЁ, (В), [= С(в х [4,5)) ‚ то из (9) легко следует, что 


+ 
ЕЕС (В ) ‚ откуда, в частности, получаем, что Е(0)=0. 
Следовательно, на основании теоремы Дини можем утверждать, 


что сходимость последовательности ТЕ, (х)} о В каждом компакте 
и= 


+ 
из К’ равномерна. 


Итак, доказана следующая 
Теорема 2. При условиях А)-С) и (14) уравнение (9) имеет неот- 
рицательное (нетривиальное), непрерывное и существенно ограничен- 


+ 
ное на ®` решение, причем, имеет место неравенство (50). 


3.2. Разрешимость основного уравнения (1) 
Сперва заметим, что функция Г(х), задаваемая формулой (10), 


является нетривиальным непрерывным и ограниченным на К реше- 
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нием уравнения (1), причем, с учетом (50) и (10) она удовлетворяет 
следующему двойному неравенству: 

—1< 1 (х) <. хЕВ. (51) 
Ниже убедимся, что при дополнительном ограничении О) на функцию 
а (х,5) уравнение (1) обладает также однопараметрическим семейст- 


вом непрерывных и ограниченных решений { я (х)} вида 


Л, (х)= Л(х+тТ), 


где Г — основной период функции &(х,5) (см. условие О)). 


пЕМЬ{0} 


Действительно, учитывая (5), из (1) будет иметь 


кола = ао даша) = 


| 
вм и ос 
 —_ 


ох, ТТ [1+пТ) ао (5) ыы 
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(хит ве т (+пт) (5) =. 


а(х+тТ ет (у) ду (5) = ГР (х+иТ) = ДР(х). 


5 


Итак, доказана следующая 
Теорема 3. При условиях А)-О) и (14) уравнение (1) обладает 


однопараметрическим семейством нетривиальных решений вида 
р, (х) — Л(х+ пТ), п=0,1,2,..., где Г(х) задается согласно фор- 
муле (10), удовлетворяет неравенству (51), а К (х) — неотрицатель- 


ное нетривиальное и существенно дграниченное решение уравнения 


(9) и удовлетворяет неравенству (50). 
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3.3. Примеры функций @(х, 5) 


Примерами функции (х,5 ) ‚ удовлетворяющих условиям А)-О), 
могут служить следующие: 


° @(х,5)= В(1+=(5) т? х), 


где 2(5 ) — непрерывная на [а4,6) функция, удовлетворяющая 
двойному неравенству: 
=(5)>0, 5 [а,Б); зир =(5)<1 
5Е[а,Б) 
* @(х,5) = (со? ха (5) (т? ха, (5), 


где ё;е С[а,Ь), причем 2, ($) > В, зир =(5) А. 
5 а,Р) 
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АННОТАЦИЯ 


В данной статье рассматривается задача нелинейного флаттера анизо- 
тропной прямоугольной пластинки, обтекаемой сверхзвуковым потоком га- 
за. Задача решается приближенно численно-аналитическим методом Галер- 
кина. Задача исследуется как в линейной, так и в нелинейной постановке. 
Произведен анализ полученных результатов с известными в литературе. Раз- 
работана расчетная программа, с помошью которой можно сделать общий 
анализ устойчивости рассматриваемой физической задачи. 

Ключевые слова: анизотропная пластинка, сверхзвуковой поток газа, 


устойчивость. 


1. Рассмотрим тонкую анизотропную прямоугольную пластинку 
постоянной толшины й. Пусть материал пластинки подчиняется обоб- 
щенному закону Гука и в каждой точке имеет три плоскости упругой 
симметрии, главные направления которых совпадают с направлениями 


ортогональных координатных линий 0©,В,)/. Координатная поверх- 





* 
Работа выполнена на кафедре Математики и математического моделирования 
РАУ врамках проекта, финансируемого Министерством образования РФ. 


42 Математическое моделирование устойчивости анизотропной прямоугольной пластинки 





ность 0,В совпадает со срединной поверхностью пластинки. Пусть 


пластинка обтекается с одной стороны сверхзвуковым потоком газа с 
невозмущенной скоростью и, направленной вдоль оси 00. Иссле- 
дуется устойчивость рассматриваемой анизотропной пластинки в по- 
токе сверхзвукового газа. 

С целью математического моделирования поставленной задачи 
принимаются следующие предположения [1-4]: 

а) Гипотеза Кирхгофа-Лява о недеформируемых нормалях; 

6) Основные предположения теории гибких пластин, считая, 

что нормальные перемещения сравнимы с толщиной пластинки; 

в) Избыточное давление газа представляется по приближенной 
формуле «прощневой теории». 

На основе принятых предположений получается следующая 


система дифференциальных уравнений движения пластинки [5—7]: 

















д“Е д“Е д“Е 
од + (в 24) то “ера 
9?» д?» д? д?» ди , о 
РА о =0, 
ЭВ до” до’ дв дооВ 
4 4 4 2 2 2 2 
и НЫ 
да до`оВ 98° до’ 98” ЭВ’ ды 
9» д?Е ® 
+7 (о,В,#) =0, 
додВ дооВ 
Здесь 
Й3 
р, о си = Ву Й, 
с с [6 
о, О Оо 
1 
Я = Ось св, 
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где и(0,В,г) — прогиб пластинки; В, — коэффициенты упругости; 
Е=Е(0,В,!) — функция напряжений. Отметим, что здесь коэффи- 
циенты анизотропии а», с„,В„ приведены согласно обозначениям ра- 
бот. Поперечная нагрузка 7 (0,В,г), возникающая вследствие возму- 


щений, складывается из сил инерции, сил демпфирования и избыточ- 


ного аэродинамического давления Др=р-р..: 


9?и ди? 
й + Др, 3 
д? Ро Е м ( ) 


где р. — плотность материала пластинки, & — коэффициент линейного 





Р=-рой 


затухания, р - время, где р — давление газа на поверхности пластин- 
ки, р.. — давление невозмущенного потока газа. 


Гидродинамическое давление В на поверхности тела вычис- 


2= 


х е-Т у. \=1 

ляется на основе формулы поршневой теории р=р.. о , 
а. 

где у, - нормальная составляющая скорости точек поверхности 


> —1 
пластинки. Разлагая это выражение в ряд Тейлора по степеням уа. и 
ограничивая членами не выше третьей степени, получим: 


2, (+!) № . &(=+1) У, - 


+ 
а 4 а 12 а 


со со со 





р=р.|1+ 


ду? ду? 


Далее, подставим вместо Уз выражение у, = —+0 — „где 


дх 


и = и(х › У) прогиб недеформированной срединной поверхности. 


Нестационарный член в правой части последней формулы соотвест- 
вует аэродинамическому демпфированию, влияние которого в боль- 


шинстве случаев оказывается незначительным, и учитывается лишь 
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> > > —1 
наиболее существенный линейный член жр.а_ д\/9#. В результате, 


формула аэродинамического давления принимает следующий вид: 





&(=+1 М? ° =(=+М? : 

р=р. По ) 0% + о 
дх 4 дх 12 дх 

Подставляя последнее в выражение (3), получим следующую 


формулу вычисления поперечной нагрузки: 


2 
ао йе [роны ера. ми 
а 





РУ ь 1 [210 
нем: [9") кем?" 
4 [210 12 [210 


где М =ИЦ/а.. -— число Маха для невозмущенного потока газа; & — по- 
казатель политропы. 

В представлении (4) из членов, учитывающих аэродинамическое 
демпфирование, сохранен лишь наиболее существенный линейный 
член. 

При исследовании вопросов устойчивости к уравнениям (1) - (2), 
где Г определяется согласно (4), присоединяются также условия на 
контуре пластинки. 

2. Рассмотрим шарнирно опертую по всему контуру ортотроп- 
ную пластинку, края которой свободно смещаются в плане. Тогда гра- 
ничные условия имеют следующий вид [7]: 

при 9=0, чо=а 








ди 9?» 
ш=0, М. = 15? 12 ЭВ? = 0, (5) 
при В=0, В=Ь 
ди 9 
0-0, М.Р, 9,570. (6) 


Приближенное решение системы (1) - (2) будем искать в виде 


И.А. Варданян 





и (о, В.Е) = У Г(е)зт а - эт В, 
1=1 


й КП 
Ве =". 
7 а и, Ь 
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(7) 


Здесь Г(Г) — подлежащие определению функции времени (. Рас- 


сматривая (7), легко заметить, что граничные условия (5) и (6) удов- 


летворяются. 


Сформулированная краевая задача детально рассмотрена в рабо- 


тах [5—7] при т=2. Известны некоторые исследования при других зна- 


чениях т [8,9], но в основном в них проведены численные расчеты и 


аналитическое исследование практически отсутствует. 


Настоящая работа посвящена численно-аналитическому исследо- 


ванию сформулированной краевой задачи методом Галеркина при лю- 


бом значении т. В частности, при т=4, подставляя (7) в (1), получаем 


следующее линейное уравнение относительно функции РЁ (0,В,#): 


4 4 4 Е, 
аа —2а,.) он о ми [со8.,0-+сози,В] ^^ + 
о 


ботов? “В; 





2:2 
А [оо5д.о +9 со5А. 4+ 9005).0-605 и,В-с0$^.0-сози-В] Д.Х, + 


+2 [с0$^. 01+ сои В] > + и [-с0$А,%+4с08^.. + 


+4 с0$).,0- со иВ-—с0$7..0- сози.В] Д.Л, + 


2-2, 
а 0$. + 256031. +25с0$1.0- соз,В-—с0$7.;0%- соз,В] ХХ + 


2.2 
+ [со$^ 0+ соз.В] 1 + 9 с0$1.0+25605/.. 0+ 


+2505... сози,В-9 с0$^..- сои >В] Д.Л, + 
+ [-с0$^.,0+9 08, 0+9 с03/.,0- сои В— с0$1..0%-сози.В] ХХ, + 


+8: [с08/0-+ сз, В] {2 + 


2.2 
[о , 0+ 49 с0$/..01-+ 49 051.0. созиВ-— 0$... сои >В] ХЛ, 


(8) 
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Решение уравнения (1), удовлетворяющего условиям (5) и (6), на- 
ходим методом неопределенных коэффициентов. В это решение вхо- 
дят неизвестные функции ХД(Г), для определения которых восполь- 
зуемся уравнением (2). Подставляя (3), (7) и найденное выражение для 


Е(0,В,г) в (2) и применяя метод Бубнова-Галеркина для определения 
безразмерных неизвестных х, = (И) /й, 1=1.2,...т, в случае орто- 


тропной пластинки получаем следующую систему нелинейных обык- 


новенных дифференциальных уравнений: 
2 


4?х. @ | 
аа (ече, иен (Л, Л)=0, (1=1,2,..т) (9) 


Так как при больших значениях т приведенная система нелиней- 





ных уравнений слишком громоздкая, здесь приводим лишь линейную 


систему для случая т-4. 








ах 4х 

И. а ве 0, 

2 

с + кух, ——Кух, =0, 

_ (10) 
4“х, 





12 
Е ем, т =0, 





и я и ве =0. 
{Е 15 И 


Таким образом, вопрос устойчивости рассматриваемой аэроупру- 
гой системы сводится к исследованию устойчивости решения системы 
(10) в зависимости от значения приведенного параметра скорости У. 


Для этой цели представим частное решение системы (10) в виде: 
—_ 7 :_1Т Л 
х. = уе”, 1=14. 
Подставляя это решение в системы (10) получим однородную систему 


линейных алгебраических уравнений относительно у,. Известно, что 


такая система имеет ненулевое решение, если ве определитель равен 
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нулю. Составляя указанный определитель, получим следующее харак- 


теристическое уравнение относительно АД: 
8 7 6 5 4 3 2 = 
ам +а № +а. № а +а. № +а. № +а № +аА+а, =0 (11) 
где: 


«=Ь а=4х, а=(у+\++1+62), а,=Х(3+31 +31 +3 +42), 


ИРУ нранзи зи ++), 
дерет КУ ранее +, 














2205 
5364? А 
+ + + + 
в = + + +15 +55 ТИ 
6 АЙ и у, 10792А?у? 
+ +; +, + + } 
> ИН И И 
53642 -; 33184 УР 
+ + + + 
О 
2..2 2:2 
мы У и У У, 2) 
м. и БАЙ 16?У т. З6А?У ри АУ’ б55364"У" 
т 225 уд 9 30625 


Невозмущенная форма пластинки будет устойчива, если действи- 
тельные части корней характеристического уравнения будут отрица- 
тельными. Следовательно, условия устойчивости, согласно теореме 
Гурвица [7], записываются в виде: 

М О) (13) 


где А, — диагональные миноры следующей матрицы Гурвица 
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аа 00 00 
И” оо 
Я. Ч: 3 оо 
Я в ве 29 2:20 
Пса. аи 34:2.10> 0 
Я’ у 2 чб 0 

9. ‚а. За: №. "О 





ба ма аа 


Система (10) решена численно при следующих исходных 


данных: Е=7.3.10°Н/м°; и = 0.34; р, =2.79. 10°кг/м” (дюралюминий), 


&=1.4;р. =1.29кг/ м° (воздух). Критическим называется то значение 


параметра скорости У =У 


сг? 


при котором нарушается хотя бы одно не- 


равенство системы (13). В случае квадратной пластинки (ф=а) для раз- 


личных значений приведенной толщины (Б/а) в Таблице 1 приведены 


значения критической скорости для двух значений т (т=2 и т=4). 
































Таблица 1. 
значение критической скорости У =У,, 
Б/а ш=2 п=4 
1/50 2.6021 3.736 
1/70 0.6773 0.9725 
1/100 0.0158 0.0222 
1/200 0.0102 0.0146 
1/300 0.002 0.0029 








Результаты вычислений показывают, что 1$) чем тоньше пластин- 


ка, тем меньше различаются результаты; 2) случай т=2 дает понижен- 


ное значение критической скорости. 
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4. Перейдем к исследованию нелинейной задачи. В этом направ- 
лении целью настоящей работы является составление расчетной прог- 
раммы-пакета в среде Мафетайса, которая для любого количества 
слагаемых в представлении Галеркина (7): 1) вычисляет функцию нап- 
ряжений, 2) получает нелинейное уравнение в частных производных, 
3) применяет метод Галеркина, 4) выдает нелинейную систему обык- 
новенных дифференциальных уравнений (9) в раскрытой форме; 5) 
вычисляет значение критической скорости, согласно критериям Гур- 
вица (13), на основе соответствующей линейной задачи. Следует отме- 
тить, что в литературе нам неизвестно существование такой програм- 
мы, а иногда ссылаются на некоторые существующие пакеты, которые 
численными методами решают задачи Коши. То есть с их помощью 
исследуется поведение амплитуды флаттерных колебаний во времени 
и все результаты, полученные при применении указанных пакетов, 
приводятся в виде графиков или таблиц. Возникает вопрос: в чем 
преимущество разработанного нами пакета? Дело в том, что мы пред- 
лагаем универсальную доступную программу, которая, в отличие от 
уже известных, дорогостоящих и недоступных для нас программ, 
имеет открытый код и имеет возможность исследовать как зависи- 
мость амплитуда-скорость, так и зависимость амплитуда-частота 
флаттерных колебаний, которые играют важную роль при исследова- 
ний задач устойчивости прямоугольных тонких тел в сверхзвуковом 
потоке газа. Корректность результатов программы неоднократно про- 
верена и подтверждена сравнением с результатами классических и 


современных авторов. 
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АННОТАЦИЯ 


В данной статье получены необходимые и достаточные условия для 
сравнения с весом одного класса многочленов двух переменных. Цель работы 
— исследовать вопрос сравнения с весом многочленов от двух переменных. 

Ключевые слова: однородный многочлен, сравнение с весом многоч- 


ленов. 


Будем пользоваться следующими стандартными обозначениями: 


ый п 
М№- множество натуральных чисел, № =№\/10}, №, пЕ М, мно- 
жество и-—мерных мультииндексов, т. е. точек © = (4,,....@,,), М, Е №, 


71=1...п, А№"- п-мерное вещественное эвклидово пространство 


точек &=(&,...&,), В =|6е Пе 20}. Для любых бЕА" и 
@Е М, обозначим &“=6“..2“, О’=рм..0%, где Е 

й 
еб 





* 
Работа выполнена в рамках тематического финансирования РАУ из средств 
МОБН РФ. 
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Характеристическим многогранником (х. м.) конечного набора 
АСВ, = [4 ЕА”",ё а 0, = Виза называется минимальный выпук- 
лый многогранник 5% < А”, содержащий множество А 210}. 


Многогранник % с А” называется полным, если 5 имеет вер- 
шину в начале координат и отличные от начала координат вершины на 
каждой оси координат. Полный многогранник Жс К, называется 
«правильным» (вполне правильным (в. п.)), если компоненты всех 
внешних (относительно Ж) нормалей (п-1)-мерных некоординат- 
ных граней неотрицательны (положительны). 


Для в. п. многогранника Ус А!" через Л(%) обозначим мно- 
жество нормалей д=(4,.../,) (п-1)-мерных некоординатных гра- 
ней ХХ, для которых ши, = 1... п} =1. 

Через В(п) обозначим множество тех в. п. многогранников 
Хс А’, для которых 


пах (у, 2) = ИА, +...+и, „УЕ Ж} = 4 (1) <1 УЕ Л(У). 


Для многогранника Хе В(п) обозначим °\°- множество вершин 


УХ и составим функцию 


(=> =. 


п 
УЕХ у 


у» 








Числа 


4, (%)= пах [4 (1), ле ^(%)} у 


р (®)= тах [И ЕИ +...+И,ИЕ 93° Е 





р, (9%) = ша, 0+уЕ 9] 
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Называются, соответсвенно, формальным, максимальным и минималь- 


ным порядками функции Й; . 


Пусть Р(б) = У 7,5“, 9. ЕС, многочлен, где сумма распрост- 


а=(Р) 
раняется по конечному набору (Р)= [ме миа 0} . Представим мно- 


гочлен Р в виде: 


Р(2)=УР(®. (1) 


1=0 





где т = плах {|0 ‚Е (Р)} — порядок многочлена Р,а Р, — однородный 
многочлен порядка ], /=0,....т. 


Точка те А" называется нулем многочлена Р порядка ((т)е М, 





если 
Р®) (т) = (2°Р)(т)=0 Уае №, |< ((т) и 
у [2 (2 #0. 
е|=Й+) 


По определению при Р(7)*0 будем считать, что ((7)=0. Для 


однородного многочлена Р обозначим 


У(Р)={те В",|1|=ЪР(т)=0}. 








В дальнейшем будем считать, что и=2 иё=(&,6.). 


Определение 1. Скажем, что многочлен О0(&) № (\еВ(2)) 
слабее многочлена Р(ё) и запишем О<”Р, если существует по- 
стоянная С > 0 такая, что 


(Е, (2))<С-Р(Ель (&)) Уе В, 


где для данного многочлена 4 и числа #> 0 
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99 
95" д5>” 


Для многоугольника Же В(2) через (1.(%),0) и (0,1, (%)) 





ЕЕ 


Е № 


9(5) 








обозначим вершины Хх ‚ Лежащие на соответствующих осях коорди- 
нат. 


Определение 2. (см. [1]) Многочлен Р(&)=»`Р‚(&) порядка т 


И 


называется гиперболическим (по Гордингу) относительно первого 


компонента, если Р, (1,0,....0)*0 и существует постоянная Се К 


такая, что Р(&)*0, 3 ЕС, (&,,...,&,)е В" при Шаё <С. 
Определение 3. ( см. [9] ) Многочлен Р(&) УР (©) порядка 


т называется 5 (5>1) гиперболическим относительно вектора 


ПЕК", если Р, (1) *0 и существует постоянная С >0 такая, что при 





1 
Ц=с[+ : 1ЕК, Р(ё+11)=0. 


Определение 4. Пусть #1(&)>1, &е А". Многочлен 


Р(2)=».Р(&) порядка т называется й гиперболическим относи- 
1-0 


тельно вектора 77е А", если Р, (1) +0 и существует постоянная С >0 
такая, что при |!|> С#(5), (ЕК, Р(6+и1)=0. 
В.Н. Маргаряном в работе [6] доказана следующая теорема. 
Теорема 1. Пусть Хе В(2), р. (Х)<1, Р, - однородный мно- 
гочлен порядка т,. Многочлен О порядка т, представленный в виде 


(1), /х слабее многочлена Р, тогда и только тогда, когда 


то 
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>, 

2) т‚-((т) (1-6, (т)) > шах {к -4, (т) (1-4, (т). к =0,...т], 

где ((т) и (,(т) порядки нулей многочленов Р, и О, 
(1=0...., т) в точке теУ(Р, },а 


р. (5%),те У Виа, 
(в) =] ие УР, =0 
в) У (р. =0 


Теорема 2. Пусть однородный многочлен порядка т Р,(&,&), 


Р,Р,-Б,Р, =0, гиперболичен относительно вектора И=(1,0), 


т 


Хе В(2), р(%)<1, О = О (5,5), где О, однородный мно- 


7=0 
гочлен порядка /, /=0,...К. Многочлен О йЙу слабее многочлена Р, 
тогда и только тогда, когда О =" Р,, где 
> и и 
Х=иЕ 2, +2 <. 
^ж(®) (У) 
йа Ро 


Доказательство. Так как +(1,0)в У то для любого 





Ме В(2) с р, (97) <! условие О<"Р, в силу теоремы 1 эквива- 
лентно следующей системе неравенств: 


к<т 
т-((т,Р,)(1- (97) > 1-6 (5,0, (1-2, (9%)),1=0,..К,те У (Р,)ПЮ @) 
т-((т,Р,)(1- 2, (9%) > 1-( (5,0, )(1- 2 (95) ,=+(0,1)е У (Р,). 
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Так как в силу определения % (%) =#(%) ир, (5) = р, (5), 
то утверждение теоремы непосредственно следует из неравенств (2) 
для 9% и % всилу теоремы 1. 


Следствие 1. Если при условиях теоремы +(1,0)е »`(Р,), то 


т 


Бо 
О=< р тогла и только — тогда, когда О—< *Р где 


хе 52, Ш = 
р. (%) 


Доказательство проводится аналогично теореме 2. 





Следствие 2. При условиях теоремы если А <т, то многочлен 


Р,+О й; гиперболичен относительно вектора 77 Е К". 


Доказательство непосредственно следует из теоремы 2 в силу 
работы [5]. 
Следствие 3. При условиях следствия 1 если К<т, то мно- 


гочлен Р+О 5=Ир.(%Ж) гиперболичен относительно вектора 


ПЕ К". 
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АННОТАЦИЯ 


Рассмотрены вопросы разработки математической модели системы уп- 
равления квадрокоптера в операторной форме с учетом динамики бесщеточ- 
ных электродвигателей постоянного тока. Показано, что введение в систему 
матричного развязывающего регулятора позволяет исключить жесткие 
взаимные связи между каналами управления по высоте полета и углам танга- 
жа, крена и рыскания, и свести проектирование системы управления к проек- 
тированию отдельных каналов на основе методов и подходов классической 
теории регулирования систем с одним входом и выходом. 

Приведен численный пример расчета системы управления квадрокоп- 
тера. 

Ключевые слова: кинематика квадрокоптера, взаимные связи, многомер- 


ная система управления, внешние возмущения, развязывающий регулятор. 
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тельской деятельности РАУ. 
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1. Введение 


Вопросам разработки математических моделей и систем управ- 
ления квадрокоптеров посвящено большое количество научно-техни- 
ческих публикаций [1-8]. При этом обычно предполагается, что бес- 
щеточные электродвигатели постоянного тока (моторы) описываются 
безынерциоными уравнениями, т.е. тяга пропеллеров пропорциональ- 
на приложенным к электродвигателям входным напряжениям. 

В настоящей статье рассмотрены вопросы разработки уточнен- 
ной математической модели квадрокоптера с учетом электромехани- 
ческой постоянной времени электродвигателей (т.е. последние 
представляются в виде апериодического звена первого порядка), а так- 
же предложен подход к проектированию системы управления квадро- 
коптера с позиций теории многосвязного регулирования. Указанный 
подход состоит во введении в систему так называемого развязываю- 
щего матричного регулятора [9-11], который компенсирует взаимные 
связи между отдельными каналами управления и сводит проектирова- 
ние взаимосвязанной системы управления квадрокоптера к проектиро- 
ванию отдельных каналов на основе методов и подходов классической 


теории регулирования систем с одним входом и выходом [13, 14]. 
2. Динамика квадрокоптера 


Схематическое изображение квадрокоптера показано на Рис. 1, 
где через О’'Х'У'7' обозначена инерциальная система координат 
(СК), относительно которой будет описываться движение центра масс 
(ЦМ) квадрокоптера, а через ОХУЙ -— жестко связанная с квадрокопте- 
ром СК с началом в ЦМ и осями, направленными вдоль главных мо- 


ментов инерции. 
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Рис. 1. Схематическое изображение квадрокоптера: 
а) системы координат, б) расположение моторов. 


Обозначим т массу квадрокоптера, х, у, 2 — его координаты от- 
носительно инерциальной СК, Г, [ ыз 1, моменты инерции относи- 


х 


тельно осей ОХ, ОУ и ОД, а 6, ф, и ц - углы тангажа, крена и рыс- 
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кания [1-4]. Далее, обозначим К,, К,, К., К», К,, К, коэффициенты 
лобового сопротивления при поступательных и вращательных движе- 
ниях квадрокоптера, „Л, — одинаковые моменты инерции ротора каж- 
дого мотора (с пропеллером), % — суммарную скорость вращения 
пропеллеров: 

0=-0-О,+О0,+0,, (1) 
где ©, — скорость 7-го пропеллера. Тогда динамические уравнения 
поступательного движения ЦМ в инерциальной СК и углового движе- 


ния квадрокоптера относительно ЦМ в связанной СК имеют следую- 


щий вид [1-5]: 














и и. (создут сои +т ту) +7. К, 

т (о =и, (с0$фзт 9 зш у — зт фсозу) + 7, —К, и (2) 
ии (созыв совр) теч, 1. 
о вой 
И: ЗЕЕ 
ча, Ч, 


где и, — подъемная сила вдоль оси вертикальной оси 0'7' ‚ < — грави- 


6? и, — управляющие моменты вокруг 


тационная постоянная, и», и 
главных осей инерции квадрокоптера, а 7,, 7,, 7, и 7,, 7,, 7, - 


внешние возмущающие силы и моменты, вызванные, например, поры- 


вами ветра. Отметим, что все переменные в уравнениях (1)—(3) яв- 
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ляются функциями времени Е, однако для краткости записи зависи- 
мость этих переменных от {Е в явном виде не указывается. 

При малых угловых скоростях квадрокоптера можно в (3) пренеб- 
речь нелинейными членами, содержащими произведения угловых ско- 
ростей, и членами в (2) и (3), зависящими от сопротивления воздуха. 
Кроме того, можно считать углы тангажа 9 и крена ф малыми, а так- 
же принять, что в условиях прямолинейного движения угол рыскания 
равен нулю (у =0) и в направлениях осей ОХ' и ОУ' в (2) приблизи- 


тельно выполняется условие и, = те . Это приводит к следующим уп- 


рощенным линейным уравнениям движения: 
2 
И 9+7 
т—-=те г 
Не ь 








а 

тт =-т8ф+У,. (4) 
422 

Е о 
а? 

7. @? = + 7ь, 
46 _ 

в аР И, +7», (5) 
Фи 

у и. НУ. 


Нелинейные уравнения движения (2) и (3) будут использованы 
для динамического моделирования системы управления квадрокопте- 
ром в среде пакета ЗипаНик. 

В современных квадрокоптерах обычно используются бесщеточ- 
ные электродвигатели (моторы) с внешним ротором [1]. Если пренеб- 


речь индуктивностью якоря электродвигателя, то тяга Т, ‚ создаваемая 
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каждым 1-м мотором, связана с входным напряжением и, дифферен- 


циальным уравнением первого порядка: 


ат 
аля —@,Т+Ки@ции,, 


(6) 
=1,2,3,4, 


где постоянные коэффициенты К м И и являются конструктивными 


параметрами, и для простоты считается, что все моторы одинаковы. 


Величина @,, в (6) при этом обратно пропорциональна электромеха- 
нической постоянной времени мотора Ту, т.е. @, =1/Т„. Отметим, 


что дифференциальные уравнения (6) фактически описывают динами- 
ку электродвигателя, т.е. определяют в первом приближении зависи- 


мость между входными напряжениями и, и угловыми скоростями вра- 
щения ОФ ,. В свою очередь, тяга Т., строго говоря, пропорциональна 
> 2 
квадрату угловой скорости ©, т.е. Т=си@,, где си =сопяЕ [1-5]. 
Однако, учитывая, что номинальная скорость вращения моторов обыч- 
но достаточно велика, то в окрестности рабочей точки квадратичная 
зависимость между Т и @„ может быть заменена приближенной 
линейной зависимостью ДТ =4„АО,, где 4» есть угловой коэффи- 


циент касательной к квадратичной зависимости в рабочей точке. Это в 


итоге приводит к уравнениям (6), где вместо АТ, используется для 


краткости обозначение Т., а коэффициент 4„ учитывается в Ку. 


3. Геометрия (кинематика) квадрокоптера 


Моторы с пропеллерами установлены на квадрокоптере таким 
образом, что задний и передний моторы (1 и 2) вращаются по часовой 


стрелке, а левый и правый (3 и 4) — против часовой стрелки (Рис. 16). 
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Каждый 1-ый мотор расположен на расстоянии Г от центра масс О и 


при вращении производит крутящий момент 7,. Создаваемая {-ым 
пропеллером тяга Т всегда направлена вдоль оси ОЙ, т.е. парал- 
лельна оси вращения роторов. Тяги Г, и вращающие моменты 7; при- 
водят к подъемной силе и, вдоль оси ОЙ и управляющим (вращаю- 


щим) моментам и,, и, и и, по осям ОХ, ОУ и ОХ, которые опреде- 


ф 
ляются следующими выражениями: 
и.=Т+Т,+Т,-+Т, 
и, = МТ -Т,), 
(7) 
и, = ИТ -Ть), 
и, ЕТ, +Т, -,-1.. 
Если считать, что пропеллеры установлены непосредственно на 


валу мотора (без редуктора), то момент 7, 1-го мотора непосредствен- 
но связан с тягой 7; соотношением 

т =К,Г, (8) 
где К» = соп5{ есть постоянный коэффициент [1-5]. 


С учетом (8) уравнения (7) можно переписать в матричной форме: 


ие т а Е 
и Е -Г 0 от 
и 0 0 ву ИЕ 
и, К, К, -К, -К, || Т, 


(9) 


откуда непосредственно получаем в матричной форме связь между 


тягами Т,Т,, Т,Т. и требуемыми управляющими моментами их, и,, 


М 
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р ле о Их 








д 
4 27 4К, 
1 1 1 и 
т” ® Зе 
и || № 
= , 10 
ь в а — 
4 22 4К, и, 
1 с 
|4 22 4К, | 


С учетом (7) и (8) уравнения (4) и (5) запишутся в виде: 








2х 1 
аг гы 7 
4? 1 
= #87, (11) 
422 1 1 
а. 
т 
420 Г 
1 (Т, 7+7 № 
а Г 
1 пе -т)+-- у, (12) 


Т, 


п, К 
р те+т- пы -т)+- -7,- 


аг 


что дает совместно с уравнениями (6) полную систему динамических 


уравнений квадрокоптера с учетом динамики (электромеханической 


постоянной времени) моторов. 


Основной особенностью системы управления квадрокоптера яв- 


ляется то, что имеется только четыре управляющих сигнала, а именно: 


напряжения и, и,, и, и и, на входах моторов, причем, тяги Г, всех 


четырех моторов направлены параллельно оси ОЙ. Подобная конфи- 
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гурация не дает возможности управлять непосредственно всеми 
шестью степенями свободы квадрокоптера (тремя линейными переме- 
щениями х, у, 2 в инерциальной СК и тремя углами крена 9, тангажа 
@ и рысканья и вокруг осей связанной СК). Поэтому в качестве четы- 
рех управляемых параметров движения квадрокоптера обычно выби- 
раются углы ©, ф, и и высота 2 [1]. В соответствии с первыми двумя 
уравнениями в системе уравнений (11), управление перемещениями х 
и у в инерциальной СК, т.е. перемещениями в горизонтальной плос- 
кости О'Х'У', при этом осуществляется косвенно, за счет изменения 
углов крена @ и тангажа ф, а управление вертикальными перемеще- 
ниями 2 вдоль оси О'7' происходит за счет изменения суммарной тя- 
ги пропеллеров и, =Т1 +Т, +Т, +Г, [см. уравнения (7), (8) и третье 
уравнение в (11)]. 

Следует указать еще на одну важную особенность системы уп- 


равления, обусловленную выбором управляемых параметров движе- 
ния квадрокоптера. Как видно из уравнений (11) и (12), на динамику 


полета оказывают влияние также внешние возмущения 7,, 7,, 7. и 
7, 7», 7, При этом только четыре из них, а именно, 7. и 7», 7, 7,, 
могут быть в некоторой степени скомпенсированы системой управле- 
ния за счет применения специальных средств. Что касается возму- 
щений 7, и У, по неуправляемым осям О'Х” и О'У' инерциальной 


СК, скомпенсировать воздействие этих возмущений за счет управляю- 


щих сигналов и, и,, и, и и, , в принципе, невозможно. 


4. Описание динамики квадрокоптера в операторной форме 
Переходя в линеаризованных уравнениях (6), (11) и (12) к 
операторной форме записи (при нулевых начальных условиях), 


получим: 
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са гевеи Я) 














(13) 
1 
у= [9+ к 
52 
1 
=== [8 +1,+1,+Т)-те+7.]. (14) 
Тед 1 
о анти 
ИЕ 7+7 (15) 
52 [45 Г, #} 
1 [К 1 
и еаа-в-тнт | 
ТЕ, (и, 
(16) 
О д. 
где 
ил, (5) =К— 9 (17) 
5 @м 


есть передаточная функция идентичных моторов. 

Отметим, что все переменные в (13)-(16) зависят от оператора 
Лапласа 5, что также в явном виде не указывается для краткости 
записи. 

Подставив (16) в (14) и (15), найдем следующие выражения, свя- 


зывающие в операторной форме переменные 2, 6, 9, и у с 


напряжениями и, и,, и, и и, на входах моторов и возмущениями 7, 


ее 


| 1 
2= 5 
т5 





1 1 
+и, +. +.) -— ны В (18) 
5 т5 
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Е 
и (5) и, (19) 
ТРЕ е. 
тр 
ф=— м (5) (м; — о т Ей (20) 
1,5 1,5 
Ку 1 
И и Е = 7 (21) 


Если ввести в рассмотрение векторы регулируемых переменных 


> ‚ Е Я 

И] = [8.8 ф, и] ‚ управляющих напряжений и =[и,и,,и.,и,| и внеш- 
<“ бека Г <“ 

них возмущений Е | ‚ а также постоянный вектор 


Е = 00 о], то уравнения (18)—(21) можно записать в следующей 


матричной форме: 








ыь \. _ 1 
Й =И,, (5) Ки ЕЕ : (22) 
где матрицы Й/, (5), К и И’, (5) имеют вид: 
Е 0 0 0 
т 
0 й 0 0 
($) = и (5) . у (23) 
0 0 — 0 
те 
К 
0 0 0 —. 
| у - 
1 1 1 1 
1-1 0 0 
В = Е (23а) 
0 0 1 - 
1 1 -1 
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Е оо 0 
т 

0 й 0 0 
И’, (5) =— ” 24 
„(5) 52 0 0 ЗВ 0 ( ) 

Т, 
ооо — 
И 








Отметим, что матрицы И, (5) и И’, (5) в (22)—(24) являются диа- 


гональными, а числовая матрица К (23а) характеризует рассмотрен- 
ные выше кинематические особенности квадрокоптера [ср. эту матри- 
цу с матрицей в (9)]. 

Таким образом, рассматриваемая система управления квадрокоп- 
тера относится к многомерным (многосвязным) системам управления, 
где взаимные связи между отдельными каналами характеризуются 
матрицей жестких (т.е. безынерционных) взаимных связей К (23а). 

Следует указать также еще на одну существенную особенность 


системы управления квадрокоптера. Подъемная сила и, вдоль верти- 
кальной оси ОЙ в (7) зависит от суммы всех тяг Т,. Однако, как выте- 


кает из вида уравнений (7), (14) и (15), если использовать одинаковые 


приращения АТ всех тяг Т › ТО это приведет только к изменению и, 


и не окажет никакого влияния на изменение вращающих моментов 


вокруг осей ОХ, ОУ и ОА ‚ так как моменты и,, и, и и, в (7) зависят 


9 
попарно от разности соответствующих тяг. С другой стороны, если 
при управлении поворотами вокруг осей ОХ, ОУ и ОХ использовать 
одинаковые по величине приращения тяг +ЛАТ, соответствующих мо- 
торов 1,2 и 3, 4 (Рис. 16), то это не приведет к изменению суммарной 


подъемной силы и,, те. к изменению высоты полета квадрокоптера. 
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Фактически, применение указанного подхода позволяет развязать 
управление квадрокоптера по всем четырем каналам управления. 

К сказанному можно подойти с позиций общей теории многосвяз- 
ного регулирования [9-12]. Допустим, ставится задача выбора такого 


матричного регулятора Й’,.,(5) замкнутой системы управления квад- 


рокоптером, который обеспечит развязку всех четырех каналов взи- 


мосвязанной системы, описываемой выражением (22). Если обозна- 
= ув 
чить через 7. =[2,,0.,ф..и.| вектор задающих (входных) сигналов 


системы управления, то матричное уравнение разомкнутой системы 
управления квадрокоптера с регулятором Й’,.,(5) примет вид: 
=. че 
Й =ИЙо($)Е+ Ру) - РЕ (25) 
где ==Й.-й есть вектор ошибки замкнутой системы, а 
И’, (5) =, (ВИ, (5) (26) 


— передаточная матрица разомкнутой системы. 


Из (26) очевидно, что если выбрать И’,.,(5) из условия 


Й’,.. (5) = В ‘а И’, ($) й (27) 
где матрица В" имеет вид: 
Гы г _ 
ЕГО 1 
4 2 4 
1 1 1 
ее а 0 Е 
В = 2 : (28) 
5.” 0 В 
4 2 4 
1 0 с 
4 2 4 
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а 481’. .(5)} есть диагональная матрица регуляторов у». (5), 
и, ($), \,(5), м,($) в отдельных каналах, то передаточная матрица 


И’, (5) (26) в развернутом виде примет вид: 





— и’. (5) 0 0 0 
т 

0 — (5) 0 0 

1 х 
И’, (5) = — мм ($) 1 ‚ (29) 
ь 0 0 — и, ($) 0 
1 
У 
К, 
0 0 0 т и’, ($) 





т.е. станет диагональной. Физически это означает, что отдельные кана- 
лы системы оказываются независимыми. В теории многосвязного ре- 
гулирования и управления регуляторы типа (27) называются «развязы- 
вающими» [9-12]. 


Выходной сигнал замкнутой системы 77 при этом примет вид: 


Я = Ф(5)Я, +Ф, (5), (5)7- ФР, (30) 


где 


Ф()=И +) Г (5), 


(31) 
Ф.()=И+и, Г 


есть передаточные матрицы замкнутой системы по выходу и ошибке, 
а Г — единичная матрица размера 4х4. 
Подставив выражения для И’, (5) (29) и И’, (5) (24) в (30), полу- 


чим [поскольку все матрицы в (30) являются диагональными] систему 


следующих четырех уравнений: 
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р (5)иу (5) т 
= - за (32 
Гиды] [ивы ни) ы 5? +; (им | о 
И о »› (3) 
1,5° + [ми (5)\, (5) 15° + [м (5)\, (5) 

_ [Мм 5)\(5) 1 
ре 1,57 + (5), (5) " ы 1,5 + Гм (5), (5) 7 (94) 
[ми (5), (5) м 1 5) 


+, (5), (5) о 1.57 + [му (5), (5) г 
первое из которых характеризует линейное движение квадрокоптера 
вдоль вертикальной оси, а остальные три описывают динамику систе- 
мы управления по углам тангажа, крена и рыскания. 

Устойчивость замкнутой системы управления квадрокоптера с 


развязывающим регулятором Й’,..(5) (27) при этом определяется 


распределением корней характеристического уравнения [9-12] 

(+ (5)]=0, (36) 
которое с учетом (32)—(35) распадается на следующие четыре уравне- 
НИЯ: 


тб? + и. (5) (5) = 0, 


1.57 + [м (5) м» (5) = 0, 
(7) 
1,5 + [ми ($)и/, (5) =0, 


1.5 + [м (5), (5) =0. 

Матричная структурная схема замкнутой системы управления 
квадрокоптера показана на Рис. 2, а на Рис. 3 даны структуры системы 
управления по высоте, а также по каналу тангажа (структуры систем 
управления по каналам крена и рыскания полностью аналогичны пос- 
ледней). Укажем, что матричная структурная схема на Рис. 2 построе- 


на на основе исходных уравнений динамики системы и порядка вклю- 
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чения физических элементов в отдельных каналах системы, за исклю- 
чением того факта, что передаточные функции одинаковых моторов 
и’ (5), которые структурно должны находиться между матрицами К" 


и А, отнесены (для упрощения матричной структурной схемы) к 


передаточной матрице разомкнутой нескорректированной системы 
И’, (5) (23). 














Рис. 2. Матричная структурная схема системы управления 


квадрокоптера. 





6) 


Рис. 3. Структурные схемы систем управления по каналам 


высоты (а) и тангажа (6). 
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Таким образом, проектирование системы управления квадрокоп- 


тера с развязывающим регулятором Й’,.,(5) (27) сводится к проекти- 


рованию независимых отдельных каналов, которое может быть осу- 
ществлено стандартными методами классической теории регулирова- 
ния [13, 14]. 


5. Заключение 

В отличие от предложенного в [8] подхода к развязыванию от- 
дельных каналов системы управления за счет применения динамичес- 
кой обратной связи, предлагаемый в настоящей статье матричный раз- 


вязывающий регулятор И’,.,(5) (27) существенно проще в реализации, 
так как состоит из статической матричной составляющей РИ 
диагональной матрицы корректирующих звеньев 17.(5), \»(5), и,(5), 
и, (5), которые выбираются на основе методов классической теории 


регулирования [13, 14]. 


6. Пример 
В качестве числового примера расчета системы управления рас- 
смотрим квадрокоптер с техническими характеристиками, приведен- 


ными в следующей таблице: 


























Символ Описание Значение 
Я момент инерции относительно осих | 7.5.10 кг. м? 
И момент инерции относительно осиу | 7.5.10кг.м? 
я момент инерции относительно оси 7 | 1.3.10 кг. м? 
одинаковые моменты инерции ротора . 
У . 6.5.10 кг. м” 
электродвигателей 
т масса квадрокоптера 0.65 кг 
у, расстояние моторов от центра масс 0.23 м 
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Допустим, в систему управления введен матричный развязываю- 
щий регулятор Й’,..(5) (27), где м,(5), м»(5), м»($), м,($) описы- 


ваются одинаковым скалярным пропорционально-дифференцирую- 
щим (ПД) регулятором следующего вида: 
Кр5 0.45 


о ы (38) 
ТИ 0.0015 +1 


Ир (5) =ь + 

Параметры ПД регулятора (38) определены на основе стандарт- 
ной методики классической теории регулирования [13]. 

На Рис. 4 показана динамическая модель системы управления 
квадрокоптера в среде пакета Зйпийик, где движение квадрокоптера 
описывается исходными нелинейными уравнениями (2), (3). 

В качестве примера реакции развязанных каналов управления на 
входные единичные ступенчатые воздействия на Рис. 5 приведен вы- 
ходной сигнал канала управления по крену, который соответствует 
техническим требованиям (отсутствие перерегулирования и длитель- 


ность переходного процесса не более 3 с). Выходные сигналы осталь- 


ных каналов управления аналогичны графику на Рис. 5. 
















































































Рис. 4. Динамическая модель квадрокоптера 
в среде пакета 5йнийиЕ. 
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Рис. 5. Реакция канала крена на входное единичное ступенчатое 


воздействие. 
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УТТН БЕСОЧРЫХС ВЕСОГАТОВ ТАКТХС ТУТО АССООМТ 
РУМАМГС$ ОЕ ЕГЕСТЕСАГ МОТОВ$ 
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УДК 539.3 


УСЛОВИЯ ПОЯВЛЕНИЯ ПОВЕРХНОСТНЫХ ВОЛН 
ГУЛЯЕВА-БЛЮСТЕЙНА 


В.М. Белубекян, М.В. Белубекян, В.Г. Гараков 
Институт механики НАН РА 


урем@ехсие.сот; тБе’Бекуап@уайоо.сот;гагакоу@уапдех.ги 
АННОТАЦИЯ 


В данной статьр рассматривается задача распространения чисто сдви- 
говых гармонических волн в слое из пьезоэлектрического материала гексаго- 
нальной симметрии класса 6 тт. Предполагается, что на одной из плос- 
костей ограничивающих слой равны нулю касательное напряжение и элект- 
рический потенциал. На другой плоскости возможны различные варианты 
граничных условий. Установлены условия, при которых тригонометричес- 
кий характер изменения по толщине слоя амплитуды волны переходит в ги- 
перболический (локализованный). 

Ключевые слова: сдвиговые волны, пьезоэлектрик, электроупругость, 


локализация. 


Т.Л. Блюстейном [1] и Ю.В. Гуляевым [2] была установлена воз- 
можность существования чисто сдвиговых волн локализованных в ок- 
рестности свободного от нагружения полупространства гексагональ- 
ной симметрии класса 6 тт. Исследование распространения волн в 
плоском волноводе из того же материала показывает, что могут быть 
волны, амплитуда которых меняется по толщине слоя в виде тригоно- 
метрических функций (в частности, без учета пьезоэффекта) так и в 


виде гиперболических функций [4-8]. Т.е. в зависимости от отноше- 
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ния толщины слоя к длине волны возможен переход решения задачи 
от тригонометрических функций к гиперболическим. Появление ги- 
перболических функций указывает на локализацию волны в окрест- 
ности плоскости ограничивающей волновод, что в пределе для корот- 
ких волн дает повехностную волну Гуляева-Блюстейна. 

В настоящей статье получены простые формулы, определяю- 
щие границы перехода решений в виде тригонометрических функций 
к гиперболическим (к локализованным) для различных граничных ус- 
ловий на плоскостях ограничивающих волновод в зависимости от 


коэффициента электромеханической связи. 


1. В прямоугольной декартовой системе координат (х,у,2) 


плоский волновод занимает область —©° < х <<, ОЗ у<р, < <2<®. 


Предполагается, что материал слоя является пьезоэлектриком гексаго- 
нальной симметрии класса 6 тт. Рассматриваются уравнения элект- 


роупругости для чисто сдвиговых волн [3]: 
САУ М0 (1.1) 


где 





с = 74 (1+), Д=-А, уф (1.2) 
р = 


Са 
В (1.1) (1.2) * — компонента упругого перемещения в направле- 


нии координаты 2,с.„- упругий модуль сдвига, е, — коэффициент 


электроупругого взаимодействия, & — диэлектрическая проницаемость, 


к -— коэффициент электромеханической связи, ф — электрический 


потенциал, Л — двумерный оператор Лапласа. 
Материальные уравнения (функциональные связи) будут исполь- 
зованы в виде: 
ди’ дФ ди дф 


Вы. а (1.3) 
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где р, РО, - компоненты вектора индукции электрического поля. 


Решение уравнений (1.1) представляется в виде гармонических 


волн: 
и= Г (у) ехр(@!-кх) ‚ И=З(у)ехр(@!-кх). = (1.4) 
Подстановка (1.4) в (1.1) дает общее решение для функции 2 (у) 

# (у) = Сяйку + Есйку (1.5) 


и обыкновенное дифференциальное уравнение: 
2 


п= У 
кс?" 


Уравнение (1.6) в случае 7 >1 имеет общее решение в виде три- 


-к (1-1) =0, (1.6) 


гонометрических функций, а в случае 0<7<1 - в виде гиперболичес- 


ких функций. Требуется определить при каких граничных условиях на 


плоскостях у= 0, в волноводе и для каких отношений ТОЛЩИНЫ ВОЛ- 


новода к длине волны имеет место тот или иной вид решения. Следую 
определению волны Гуляева-Блюстейна, решение в виде гиперболи- 
ческих функций дает волны локализованные в окрестности границы 
волновода. 


Для начала решение уравнения (1.6) запишем в виде 
= АзйкЛ-пу-+ Всйк ПУ (1.7) 
и потребуем, чтобы решения уравнений (1.1) удовлетворили граничным 
условиям равенства нулю касательного напряжения и электрического по- 
тенциала на плоскости у = 0 (как и в задаче Гуляева-Блюстейна). 
0, =0, Фф=0 при у=0 (1.8) 
С учетом (1.2), (1.4), (1.5) общее решение для электрического 


потенциала будет: 
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р сылку + Рейку+[лзвкупу веку) кри к (1.9) 
[2 


Определяя 0..из (1.3), (1.4), (1.7) и требуя, чтобы вместе с (1.9) 


были удовлетворены граничные условия (1.8), получим: 
с=-9“ (1+ у) Л-ПА, Е =-9 В (1.10) 
е; = 


Оставшиеся произвольные постоянные А,В должны быть опре- 
делены из граничных условий при у=БЬ. 

2. С целью установления условий появления локализованных 
волн в окрестности свободной от напряжения границы у=0 рас- 
сматриваются различные варианты граничных условий при у=Ь. 
Пусть заданы условия закрепления края и равенства нулю 
электрического потенциала 

и\= 0, ф=Оприу=Ь (2.1) 

Нетрудно проверить,что согласно (1.2) выражению для жи и 
из (1.4), условия (2.1) приводятся к условиям: 

Х=0, =0 приу=Ь (2.2) 

Подстановка (1.5), (1.7) в граничные условия (2.2) приводит к од- 
нородной алгебраической системе уравнений относительно произ- 


вольных постоянных А, В: 
Азйк/1-76 +Всейк\/1-ПЬ = 0 (2.3) 
44 (1+ 7) Л-ИАзйкЬ+ 95 ВсйкЬ = 0 
е; = 
Равенство нулю детерминанта системы (2.3) дает дисперсионное 


уравнение, определяющее безразмерный параметр фазовой скорости 


ВОЛНЫ 7] 


уйкЛ-тЬ= (1+) Л-пкЬ (2.4) 
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Из уравнения (2.4) для большинства кф (коротких волн) в при- 
ближении 

Ик -ПЬ =, ИкЬ=1 (2.5) 

получается значение параметра фазовой скорости задачи Гуляева- 


Блюстейна в случае полупространства [1,2] 
2. 


п=1- Я (2.6) 
(1+1) 


С другой стороны, уравнение (2.4) имеет корень 7 =1. Разделяя 
уравнение (2.4) и переходя к пределу (таким образом исключая корень 
П-Ь 

ДКЬ = (1+) йкЬ (2.7) 


Если обозначить корень уравнения (2.7) относительно кБ через 


(кЪ)., то гиперболическая форма решения, т.е. решение удовлетво- 
ряющее условию 

0<7<1 (2.8) 
будет для 

кф > (кБ). (2.9) 


Численные расчеты показывают, что решение уравнения (2.7), 
полученное для достаточно коротких волн (*Ь), >1, так что можно 
принять Шхкф=1 [8]. В этом приближении получается очень простая 


формула для условия существования локализованной поверхностной 


ВОЛНЫ 


(кв) >19 (2.10) 
я 


Второй вариант граничных условий на плоскости у=Ь предпо- 


лагает равенство нулю упругого перемещения и нормальной компо- 


ненты вектора индукции электрического поля: 
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и=0, О, =0 приу=Ь (2.11) 

Согласно выражению для ШО, из (1.3) и выражению для у из 
(1.2), следует ди /ду=0. Поэтому из (2.11) получается: 

Л=0, &'=0 приу=Ь (2.12) 

Удовлетворение граничным условиям (2.12) аналогично преды- 


дущему случаю, приводит к системе однородных алгебраических 


уравнений относительно произвольных постояных А,В. Равенство 


нулю детерминанта этой системы дает искомое дисперсионное урав- 


нение 


дк -пыйкь = (1+) Л-П (2.13) 

Из уравнения (2.13) также, как и для уравнения (2.4), для боль- 

ших значений кф (для коротких волн) получается решение Гуляева- 

Блюстейна. Условие существования локализованных волн (хф), полу- 
чается (при 7—1) из решения уравнения 

ДкЫйкь =1+Х (2.14) 


которое в приближении Шхкф =1 совпадают с условием (2.10). 


3. Нетрудно показать, что третий вариант граничных условий 


0; =0, О, =Оприу=Ь (3.1) 
приводится к условиям 

И, И биреВ (3.2) 

ду ду 


Требование, чтобы решение (1.4) с учетом (1.5) и (1.6) удовлет- 
воряло граничным условиям (3.2), приводит к системе алгебраических 
(однорородных) уравнений относительно произвольных постоянных 


А,В. Равенство нулю детерминанта этой системы дает следующее 


дисперсионное уравнение 


(1+7) Л-пйк.Л-йкЬ = дйкЬ 3) 
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Уравнение (3.3) в коротковолновом приближении (кф >>1) опреде- 


ляет фазовую скорость поверхностной волны Гуляева-Блюстейна 
(2.6). Уравнение (3.3) всегда имеет решение, удовлетворяющее усло- 
вию (2.8) — условию существования локализованной в окрестности 


края у=0 волны. В этом нетрудно убедиться, если рассмотреть дру- 
гой предельный случай — приближения длинных волн (кф << 1). В этом 


случае из уравнения (3.3) для фазовой скорости получается: 
ИЯ 
п=1--^— (3.4) 
ТН 

что существенно меньше фазовой скорости волны Гуляева-Блюстейна 
(с учетом 7<1). 

Наконец рассматривается четвертый вариант грничных условий, 
когда на плоскости у=Ь удовлетворяются те же граничные условия 
(1.8), что и на плоскости у=0. 


В этом случае задача распространения гармонических волн в 
слое разбивается на две задачи: 1) антисимметричные по толщине от- 
носительно перемещения у и симметричные относительно электри- 
ческого потенциала ф, 2) симметричные относительно у’ и антисим- 
метричные относительно Фф. Для первой задачи в срединной плос- 
кости слоя должны иметь место условия 

и=0, дф/ду=0Опри у=0,56 (3.5) 

Задача решения уравнений (1.1) с граничными условиями (1.8), 

(3.5) приводится, аналогично предыдущим задачам к решению дис- 


персионного уравнения: 
и ®у=я ие = хи). 
(3.6) 


Уравнение (3.5), как и уравнения (2.4) и (2.13) имеет решения 
как при 7>1, так и при 0<7<! в зависисмости от относительной 


длины волны, кф и коэффициента электромеханической связи 7. 
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Вторая задача (симметричная относительно м’) имеет место при 


граничных условиях: 


аи ф=0 приу= 0, 56 (3.7) 
ду 


В этом случае дисперсионное уравнение задачи получается в 


виде: 


д+(1+х) Л-п5йкЛ-П Бан = дейкЛ-П ысйкУ (3.8) 


Уравнение (3.8), аналогично уравнению (3.3), всегда имеет ре- 


шение в виде тригонометрических функций (7>1). В частности, в 


2 
длинноволновом приближении -(кф) <<1, при сохранении двух чле- 


нов разложения по аргументу получается: 


2(1+27 
й= 22) Ы (3.9) 

РЕЗ 
Следовательно, для задачи слоя с одинаковыми условиями ра- 
венств нулю касательно напряжения и электрического потенциала на 


обеих плоскостях у=0 и у=Ь возможны одновременно решения как 
локализованные (0<17<1), так и нелокализованные (7>1). Задача 


имеет двойную систему собственных функций. 


Заключение. Приведены решения задач распространения элект- 
роупругих волн в слое, когда на одной из границ равны нулю напря- 
жение и электрический потенциал, а на другой заданы различные ва- 
рианты граничных условий. Установлены появления локализованных 
волн для различных граничных условий в зависимости от относитель- 


ной длины волн и коэффициента электромеханической связи. 
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АННОТАЦИЯ 


Этанольный экстракт Ррейсит регопиит Г., стандартизованный по 
гиперицину, обладает высокой антирадикальной активностью и может быть 
использован для осуществления синтеза и конгрегации металлических нано- 
частиц, в частности, ЕеМР$. Нами синтезированные наночастицы обладают 
бактерицидным свойством в отношении Езйенста сой К-12, ингибирующая 
зона которых намного раз превышает зону ампицилина. 

Ключевые слова: наночастицы железа, экстракт Шрейсит ре огайит, 


Езйетста сой, антибактериальная активность. 


Введение 
С возникновением резистентности патогенных штаммов бакте- 
рий к антибиотикам в медицине назревает необходимость содания но- 


вых систем дезинфекции. В связи с этим появляется интерес к нано- 
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довательской деятельности РАУ на 2016г. 
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частицам, которые обладают асептическими свойствами. Поиск и раз- 
работка новых подходов синтеза наночастиц с широким спектром те- 
рапевтического действия из разных источников, в том числе и расти- 
тельного, с использованием новейших технологий, имеет очень боль- 
шое значение для практической медицины. Преимущества нанотехно- 
логий растений над традиционными химическими методами синтеза 
наночастиц заключаются в том, что они не используют токсичные хи- 
микаты во время синтеза, а также легко синтезируются в больших 
масштабах [1; 2]. 

В ходе обширных исследований химического синтеза нано- 
частиц было выяснено, что эти частицы также обладают несколькими 
потенциальными опасностями, включая канцерогенность, генотоксич- 
ность, цитотоксичность и общую токсичность. Размеры инженерных 
наноматериалов — наночастиц настолько малы, что могут пройти 
сквозь кожу, легкие и кишечный тракт с неизвестными последствиями 
для здоровья человека [3]. Выбор растений и отбор экологически 
чистых нетоксичных восстанавливающих и стабилизирующих агентов 
растворителя являются наиболее важными вопросами, которые долж- 
ны быть рассмотрены в зеленом синтезе наночастиц [4; 5]. 

ретсит регрюогаит (зверобой продырявленный) — это одно из 
наиболее изученных лекарственных растений во всем мире, в листьях, 
цветках и плодах которого синтезируются и накапливаются биологи- 
чески активные вторичные метаболиты [6]. Фитофармацевтические 
препараты на основе стандартизированных экстрактов, полученных из 
цветущих верхушек этого растения, были эффективны против легкой 
и умеренной депрессии [7; 8]. Н. рейогиит и некоторые другие виды 
рода важны и используются в пище как корма, топливо, красители [9; 
10], как лекарственные препараты, используемые при ранах, ожогах, 
порезах, спазме желудка, бессоннице и мышечных болях [11; 12]. Ак- 


тивные компоненты Н. ре’/огаиит, такие, как: гиперфорин, адгиперфо- 
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рин, гиперицин, кверцетин обладают антидепрессантной, противовос- 
полительной, антибактериальной, противоопухолевой, спазмолитичес- 
кой, гипотензивной, ранозаживляющей активностью [13—15]. Фарма- 
цевтическая важность этого растения быстро возросла в последние го- 
ды, особенно в связи с обнаружением вышеуказанных свойств, а так- 


же выявлением фотодинамической активности гиперицина. 


Материалы и методы 

Сбор и хранение листьев, цветков и плодов Н.регопиит, расту- 
щих в Тавушском районе Армении, проводили в сентябре 2015 года. 
Высушенные теплым воздухом (55+1°С) сухие листья, цветки и полды 
Н. реогийит (З3мг/мл) после механической гомогенизации экстраги- 
ровали в 96% этил спирте и в течение 15 минут подвергали воздейст- 
вию ультразвуком (75Вт) (Огазотис Нотозеплтег, Зот1с-150\/, МКС). 
После 24 часовой инкубации на качалке (60—70 об/мин.) экстракты 
центрифугировали 15 минут при 3000 об/мин. на центрифуге (Зоцап 
08412) при 3000 об/мин. Экстракты были стандартизованы по содер- 
жанию гиперицина (35 мкМ). 

Антирадикальную активность (АРА) тестируемых образцов оп- 
ределяли по методам тушения стабильного радикала 2,2-дифенил-пик- 
рил гидразила (ДФПГ) (ЕШКа, Сегтапу), измеряемому поглощением 
при 517 нм [16]. 

Спектры поглощения экстрактов зверобоя и полученных нано- 
частиц были исследованы с помощью абсорбционной спектроскопии 
при длине волны 200-800 нм на сканирующем спектрофотометре 
(ЦУ/У1 ЛЕМ\МАУ 6405, ВБ). 

Содержание гиперицина в экстрактах определяли абсорбционной 
спектроскопией, принимая коэффициент молярной экстинкции при 
590 нм. равным 43500. Суммарное содержание гиперицина в высу- 


шенной ткани в процентах (Х) определяли по уравнению: 
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Х-=Ах 125 /А*х п, где: А — абсорбция образца при 590 нм, А!® 
— коэффициент абсорбции, равный 870 М/см, т - навеска сухого препа- 
рата в граммах, Х — содержание гиперицина в сухом препарате [15]. 

Для получения наночастиц железа (ЕеМР$) растительный экст- 
ракт смешивали с раствором ЕеС15*6 Н›О (1-3 мМ). В процессе конг- 
регации переход Ее’, Ее? до Ее? подтверждалось изменением цвета 
раствора. Наночастицы визуализировали с помощью светового мик- 
роскопа (АтЗсоре, США). 

Исследование антибактериальной активности проводили диск- 
диффузионным методом на агаре [17]. Индикатором антибактериаль- 
ной активности являлась зона ингибирования, формирующаяся после 
24—48 часов инкубации бактерий [15]. Во все чашки вводили 0.1мл 
взвеси культуры дикого типа Е.сой К-12, природно-лизогенного штам- 
ма (штамм из коллекции кафедры Микробиологии ЕГУ) [Уаг4ареуап 
её а|. 2012]. В качестве контроля использовали стандартизованный по 
содержанию гиперицина (35 мкМ) этанольный экстракт зверобоя, 
раствор ЕеС1-*6Н>О, минимальная ингибирующая концентрация для 
положительного контроля — ампицилина составила 5 мкг/мл, а для «зе- 
леного» ЕеМР$ — 60мкг/мл. 

Статистический анализ проводили на основе комплексного при- 
менения стандартных статистических методов: вычисления средних 
значений, стандартных отклонений, стандартных средних ошибок. В 
таблицах, на графиках и диаграммах приведены средние арифметичес- 


кие и их стандартные ошибки (п = 8-12). 


Результаты и обсуждение 

Синтез наночастиц с использованием растений очень рентабелен 
И поэтому может быть использован как экономичная и ценная альтер- 
натива для крупномасштабного производства. Экстракты различных 


растений используются как восстанавливающие и стабилизирующие 
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агенты для биосинтеза наночастиц. В связи с этим проводили предва- 
рительный скрининг этанольных экстрактов растений Н. ре’огаит, 
О. БазШсит, Г. по 5 по АРА (Рис. 1). Как видно из полученных ре- 
зультатов, высокой АРА активностью обладает этанольный экстракт 


зверобоя, который и был выбран для конгрегации зеленых наночастиц. 


100 
= Еф-АЬХ 

зо ^ Н.регГогайиа . 
о ОБазШсиа 

я 60 

= 

ы 

5 40 й 





о 9.4 0.8 1.2 1.6 2 
Концентрация (мг/мл) 


Рис. 1. Антирадикальная активность этанольных 


экстрактов разных растений. 


Проводили реакцию синтеза наночастиц, смешивая стандарти- 
зованный экстракт зверобоя с ЕеС1з*6 Н2О. Визуализировали получен- 


ные ЕеМР$ под микроскопом (АтзЗсоре, США) (Рис. 2). 
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Рис. 2. Визуализация агрегатов наночастиц (увеличение 200 раз). 


Как показали наши результаты, экстракт зверобоя, раствор соли 
Ее в темновых условиях не обладают выраженной цитотоксичностью в 
отношении Р.сой К-12 (Рис. 3). Однако нами синтезированные нано- 
частицы обладают выраженным бактериоцидным свойством, действуя 
на выживаемость РЁ. сой. Зона полного ингибирования ЕеМРз равна 28 
—32 мм (учитывая кор агара 9 мм) и с площадью - 44,1-59,2 см” (без 
кора) (Табл. 1, Рис. 3). Сравнивая полученные нами результаты с 
действием ампициллина на выживаемость Ё. сой, можно сказать, что 
несмотря на то, что МИК ЕеМР$ почти в 10 раз больше, но ингибирую- 


щая зона превышает ее почти в 1,5 раза. 
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Таблица 1. 


Антибактериальная активность стандартного и тестируемого 
препарата против бактериального тестового организма (Е. сой). 


























Езйенсша сой 
Препарат Концентрация Ингибирующая 
(мкг/мл) зона (мм) 
5 14 
Ампициллин [18] 50 16 
250 20 
е 60 9 
ЕеМР$ (зеленый) 180 32 

















Рис. 3. Образование зоны ингибирования после инъекции 


«зеленых» наночастиц железа. 
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Наночастицы представляют интерес благодаря уникальным 
свойствам, которые могут быть включены в антимикробные системы, 
биосенсорные, сверхпроводящие и криогенные материалы, компози- 
ционные волокна, косметические продукты [19; 20]. Комплексные ис- 
следования антибактериального действия наночастиц появились с 
2004 года и выросли по экспоненте. Нанометаллические системы об- 
ладают рядом преимуществ, благодаря которым их используют в ка- 
честве противомикробных агентов. Они обладают очень высокой ак- 
тивностью против широкого круга микробов и паразитов даже при ис- 
пользовании низких доз. Полученные нами наночастицы можно ис- 


пользовать в качестве бактерицидных, асептических агентов. 
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АНТИРАДИКАЛЬНЫЕ СВОЙСТВА И ДЕЙСТВИЕ 
НА ЭРИТРОЦИТЫ ЧЕЛОВЕКА ЭКСТРАКТОВ 
ОСШМОМ $ АМСТОМ И ОСММОМ АКАВАТОМ` 


Ш.А. Казарян 
Российско-Армянский университет 
зризратйк.Кагагуап@гаи.ат 
АННОТАЦИЯ 


Были исследованы общее содержание флавоноидов (ОСФ), антиради- 
кальная активность (АРА) и гемолитические свойства водно-этанольных 
экстрактов листьев Осйпит 5апсит и Отит агагйит, прорастающих на 
территории Котайкского марза Армении. Выявлены высокие ОСФ и АРА 
этанольных экстрактов, с максимальным ОСФ и АРА у 50% этанольного 
экстракта О.агагиит. Обнаружено, что экстракты обладают сильной гемоли- 
тической активностью как при темновой инкубации, так и при облучении. 
Наибольшая цитотоксическая активность наблюдается при воздействии на 
эритроциты человека 96% этанольным экстрактом О.5апсит. 

Ключевые слова. О.бапсит, О.агагйит, общее содержание флаво- 


ноидов, антирадикальная активность, цитотоксичность, гемолиз. 


Введение 
Согласно данным Всемирной Организации Здравоохранения, бо- 
лее 80% всего населения мира для поддержания своего здоровья ис- 


пользует натуральные продукты, в том числе и БАДы [1]. 


* а 

Исследование выполнено в Российско-Армянском университете (РАУ) за счет 
средств, выделенных в рамках субсидии МОН РФ на финансирование научно-иссле- 
довательской деятельности РАУ на 2016г. 
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Род Осйпит (базилик) включает более 150 видов и считается 
одним из крупнейших родов семейства Гапйасеае. Все больше воз- 
растает интерес к изучению свойств растений рода Осйпит в связи с 
широким распространением, частым использованием в национальной 
кухне и традиционной медицине разных народов и содержанием боль- 
шого спектра вторичных метаболитов [2]. В народной медицине его 
используют в качестве средства от головной боли, кашля, заболеваний 
желудочно-кишечного тракта, почек и т.д. [3]. Исследованы антимик- 
робные свойства некоторых видов базилика. Доказана его эффектив- 
ность при лечении некоторых форм вирусных заболеваний. Примеча- 
тельно, что при этом смягчается проявление побочных эффектов, 
часто проявляющихся при применении синтетических антимикробных 
препаратов [4]. Кроме этого, изучены и противовоспалительные, про- 
тивоопухолевые, антиоксидантные и противогрибковые свойства не- 
которых видов Осйпит [5]. Осйтит проявляет ингибирующую актив- 
ность в отношении ВИЧ-1 обратной транскриптазы и способен подав- 
лять агрегацию тромбоцитов, индуцированную коллагеном и АДФ [6]. 
Эксперименты на крысах показали, что экстракт листьев О.Ба5уШсит 
увеличивает активность глутатион-трансферазы, которая в определен- 
ной степени способна защитить клетки органов ЖКТ от химического 
карценогенеза [7]. 

Главными вторичными метаболитами, обнаруженными в расте- 
ниях рода Осйпит, являются эфирные масла, фенольные компоненты 
и антоцианины [8]. Основными компонентами эфирного масла базили- 
ка являются эвгенол, изоэвгенол, хавикол, эстрагол и линалоол. Эфир- 
ные масла базилика в м уЙго-исследованиях проявляют сильные ан- 
тиоксидантные, противоопухолевые, антимикробные и противовирус- 
ные свойства [5]. В базилике обнаружены также розмариновая и ко- 
фейная кислоты, обладающие сильными антибактериальными и про- 


тивовоспалительными свойствами [9]. Флавоноиды и другие феноль- 
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ные компоненты базилика обладают широким спектром свойств, а 
такой класс флавоноидов, как антоцианины, обеспечивают характер- 
ный цвет пурпурных видов данного рода [10]. Помимо этого, антоциан 
цианидин-3-гликозид, встречающийся в О.агагиит, подавляет @2/М 
фазу клеточного цикла, вызывая апоптоз [11]. 

В связи с этим, следует учитывать и тот факт, что разные виды ро- 
да Осйпит могут обладать различными свойствами и, несмотря на зна- 
чительный перечень полезных биологических свойств экстрактов 
Остит, возможно наличие и цитотоксических свойств. Наиболее ин- 
формативной моделью исследования в данном случае является изуче- 
ние цитотоксического действия экстрактов на эритроциты человека, в 
связи с чем, в данной работе было проведено сравнение антирадикаль- 
ных и гемолитических свойств экстрактов двух разных видов рода 
Остит — О.бапаит и О.агагйит, прорастающих на Армянском на- 


горье. 


Материалы и методы 


Растения были собраны в июне 2016 года в Армении (Котайкс- 
кий марз). Листья растений после предварительной промывки, стери- 
лизовали в 1% растворе гипохлорида натрия и высушивали до 10% 
влажности. 

Для экстрагирования использовали сухие измельченные листья в 


соотношении 1:30 мл. Сухие листья после механической гомогениза- 





ции подвергали воздействию ультразвуком при 75Вт и времени воз- 
действия 15 мин. (О№газоп1с Нотозеттег, Зотс-150\/, МВС). После 
24-часовой инкубации на качалке (60—70 об/мин) экстракты центрифу- 


гировали 15 минут при 3000 об/мин на центрифуге юпап @В412 [12]. 
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Общее содержание флавоноидов (ОСФ) определяли за счет 
способности флавоноидов хелатировать атомы АГ", при 430 нм, по 
методу [13]. 

Антирадикальную активность (АРА) определяли методом туше- 
ния свободного радикала ОРРН (2,2-дифенил-1-пикрилгидразил) при 
517 нм на спектрофотометре ЛЕМ\УАУ 6405 ЦУ/У!15 Зреспарвоютен 
(Великобритания) [14]. 

Для исследования гемолитической активности экстрактов значе- 
ния их рН доводили до 7.2. Кровь брали у 15 практически здоровых 
доноров и центрифугировали со скоростью 6000 об/мин. в течение 10 
мин. при 5°С, далее осажденные эритроциты дважды отмывали цент- 
рифугированием при тех же условиях в физиологическом растворе. 
Осажденные эритроциты ресуспендировали в фосфатном буфере (рН 
7.2). Оптическую плотность (ОП) суспензии эритроцтов доводили до 
0,7 оптических единиц, где число эритроцитов равно приблизительно 
6.2*10° клеток/ мл. Цитотоксическое действие растительных экстрак- 
тов на эритроциты были исследованы как непосредстевенно после до- 
бавления экстрактов к суспензии эритроцитов, так и после их совмест- 
ной инкубации в темновых условиях при 25° С в течение 0-24 час. Для 
выявления фотодинамических свойств всех экстрактов суспензию 
эритроцитов совместно с экстрактами облучали видимым светом лампы 
накаливания («Этюд», РФ) с мощностью излучения на поверхности об- 
разца 30 мВт/см” при 25°С. Резистентность (устойчивость) эритроцитов 
оценивали по изменению оптического поглощения суспензии эритроци- 
тов при 680 нм в течение 15-30 мин, снимая показания через каждые 30 
сек. на спектрофотометре (ЦУ/У15 6405 “Депуаеу”, ОК) [15]. 

Биологическая повторность опытов 4—6 кратная при проведении 
2-3 серий опытов в каждом. Статистический анализ проводили с по- 


мощью программы М1сгозой ОЁйсе Ехсе]. 
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Результаты и обсуждение 


Данные общего содержания флавоноидов водно-этанольных 
экстрактов листьев О.бапсит и О.атагщит приведены на Рис. 1. Ана- 
лиз общего содержания флавоноидов показал, что наибольшее их ко- 
личество содержится в 50% этанольном экстракте О.агагаит 
(36,34=0.8 мкг/г), а наименьшее — в 96% этанольном экстракте 
О.агагаит (2.45-0.044 мкг/г). Однако следует отметить общую тен- 
денцию: как среди экстрактов О.бапсит, так и среди экстрактов 
О.агагаит, наибольшее содержание флавоноидов наблюдается у 50% 
этанольных экстрактов, наименьшее — у 96%. За исключением эта- 
нольных экстрактов, во всех остальных случаях ОСФ выше у экстрак- 
тов О.агагщит. Этот факт может быть связан с тем, что, в отличие от 
О.5апсит, в листьях О.агатйит содержатся такие гидрофильные фла- 
воноиды, относящиеся к антоцианам, как цианидин-3-(ди-р-кумарилг- 
люкозид)-5-глюкозид и пеонидин, которые также обеспечивают пур- 
пурный цвет растения [16]. 

Максимальная АРА выявлена у 50% (0,26-0,04 мг/мл) и 25% 
(0,26+0,06 мг/мл) этанольных экстрактов О.агаищит, примечательно, 
что ОСФ в данных экстрактах было выше, чем во всех остальных. 96% 
этанольный экстракт О.агаижит хоть и обладает меныпим ОСФ, чем 
96% этанольный экстракт О.бапсит, но значение АРА у него выше 
(1,8=0,05 мг/мл), а значение АРА 96% этанольного экстракта 
О.5апсит самое низкое среди всех экстрактов (2,75+0,06 мг/мл). Зна- 


чения АРА остальных экстрактов представлены на Рис.1. 
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Рис.1. ОСФ и АРА (1С50) водно-спиртовых экстрактов О.запс ит и 
О.агагаит. 


Цитотоксическое действие экстрактов проверяли при темновой 
инкубации и облучении. В условиях темновой инкубации водные и 
25% этанольные экстракты не приводят к гемолизу, 50% и 70% эта- 
нольные экстракты О.бапсит приводят к частичному гемолизу, но да- 
же по прошествию 24 ч. инкубации гемолиз составляет 36% и 18%, 
соответственно. 


Таблица 1. 


Устойчивость эритроцитов (в процентах) при воздействии водно- 
спиртовыми экстрактами О.агагаит и О.запсит при облучении. 















































Экстракт Время облучения 
0 мин 5 10 15 20 25 30 24 ч 
мин | Мин мин мин мин | Мин 

25 % 106 104 101 97 84 82 80 10 
Е [50% | 10| 9 | 25 8 - 
Ч 70 % 102 33 6 а — 
5 | 96% 114 124 123 68 67 45 30 0 
о 

25 % 89 89 93 94 100 = 88 84 
8 150% 100 98 5 — — 
$ | 70% 90 84 40 3 _ 
$ | 96% 2 1 = в = 
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Этанольный экстракт О.5апсит сразу после добавления к взве- 
си эритроцитов начинает приводить к незначительному гемолизу, ко- 
торый к 60 мин. темновой инкубации составляет 15%, после 24 ч. ин- 
кубации наблюдается полный гемолиз. Следует отметить, что 50% 
этанольный экстракт О.агагаит, обладающий наибольшими значе- 
ниями ОСФ и АРА, уже к 60 мин. темновой инкубации приводит к 
полному гемолизу эритроцитов. 70% и 96% — этанольные экстракты 
О.агагиит приводят к гемолизу, начиная с 15 и 30 мин., соответствен- 
но, а после 24 ч. инкубации гемолиз составляет 22% и 59%. 

При облучении самой высокой цитотоксичностью обладает эта- 
нольный экстракт О.5апсит, который сразу после добавления к взвеси 
эритроцитов приводит к полному гемолизу. Высокой цитотоксич- 
ностью обладают 50% и 70% этанольные экстракты О.агагйит и 
О.5апсит, которые уже при 5—10 мин. облучении приводят к сильно- 
му гемолизу. Слабой цитотоксичностью обладают лишь водные экст- 
ракты. Остальные данные цитотоксического действия экстрактов 
О.агапиит и О.бапсит приведены в Таблице 1. 

Из всего вышеизложенного следует, что в листьях О.агагаит и 
О.5апсит содержатся как гидрофильные, так и гидрофобные флаво- 
ноиды, однако с увеличением содержания этанола в экстрагенте 
уменьшается ОСФ, и увеличивается выход некого фотодинамического 
вещества, которое усиливает гемолиз эритроцитов. При этом, содер- 
жание данного гидрофобного вещества в листьях О.5апсит намного 
больше, чем в листьях О.агагиит. В связи с неустойчивостью анто- 
цианов ко многим факторам, в том числе и свету, можно предполо- 
жить, что такими компонентами являются представители данного 
класса флавоноидов [17], однако, учитывая, что все антоцианы имеют 
синий, красный и пурпурный цвет, не характерный для О.бапит, 
можно утверждать, что гемолические свойства экстрактов обусловле- 


ны веществами других классов вторичных метаболитов. Несмотря на 
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высокие ОСФ и АРА, этанольные экстракты листьев О.агагаит и 


О.5апсит проявляют цитотоксическую активность в отношении эрит- 


роцитов человека. 
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АННОТАЦИЯ 


В современной антираковой терапии большое внимание уделяется воз- 
можности ингибирования белка теломеразы в раковых клетках посредством 
стабилизации С-квадруплексной ДНК, формирующейся на концах теломер. 
В настоящее время одной из важнейших проблем является поиск низкомоле- 
кулярных лигандов, с высокой специфичностью связывающихся с данными 
структурами для их стабилизации. 

Целью данной работы являлось проведение виртуального скрининга 
для выявления низкомолекулярных лигандов, связывающихся с теломерной 
О-ДНК, а также анализ полученных пространственно-энергетических пара- 
метров. 

Ключевые слова: теломераза, антираковая терапия, С-квадруплекс, 


виртуальный скрининг. 


Введение 
Теломеры представляют собой специализированые участки ДНК 


на концах линейных хромосом, которые обеспечивают их защиту от 
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повреждений в период клеточного деления, от действия нуклеаз, а 
также от негомологичного спаривания хромосом. Теломерные концы 
ДНК человека состоят из тандемных повторов последовательностей 
гексануклеотидов 4(ТТАССС), 5-10 кб длиной, образующих терми- 
нальный одноцепочечный 3’-оверхэнг, включающий 35—600 нуклео- 
тидных оснований [1]. Теломерная ДНК ассоциирована с белками 
шелтеринового комплекса [2], а структура и стабильность теломер 
тесно связана с раковыми заболеваниями, старением и генетической 
стабильностью [3]. 

С-квадруплексные структуры, которые формируются на дисталь- 
ном 3’- конце теломер человека и содержат гуанин-богатый оверхэнг 
с приблизительно 100—200 оснований нуклеотидов, рассматриваются 
как потенциальная мишень для новой антираковой терапии, основан- 
ной на использовании низкомолекулярных лигандов. В т уйто иссле- 
дованиях показано, что ряд низкомолекулярных лигандов, стабилизи- 
рующих теломерные С-квадруплексные структуры, ингибирует актив- 
ность теломеразы [4—6]. 

Несмотря на большое количество работ, направленных на изуче- 
ние С-квадруплексных структур, прогресс в поиске потенциальных 
низкомолекулярных лигандов, стабилизирующих С-ДНК, ограничи- 
вается лишь одним соединением, допущенным до клинических испы- 
таний [7]. Одной из причин является отсутствие структурных данных 
высокого разрешения для многих предполагаемых С-квадруплекс — 
образующих последовательностей. Другой основополагающей проб- 
лемой является динамический полиморфизм этих структур и ограни- 
ченность традиционных биофизических методов в изучении данного 
явления. Одна и та же последовательность способна формировать сра- 
зу ряд топологически различных в структурном плане С-ДНК. Малей- 
шие изменения в последовательности, вплоть до одного нуклеотида, 


могут привести к критическим структурным изменениям. Спектроско- 
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пические методы низкого разрешения, такие, как циркулярный дих- 
роизм или УФ-спетроскопия недостаточно точны для определения 
различий между структурными конформерами. Методы высокого раз- 
решения, такие, как ЯМР и кристаллография, тоже предоставляют ог- 
раниченную информацию при исследовании динамического полимор- 
физма С-ДНК [8]. 

Целью данной работы являлся поиск низкомолекулярных лиган- 
дов, связывающихся с рядом топологически различных структур тело- 


мерной С-ДНК с высокой энергией связывания. 


Материалы и методы 


В качестве мишеней были отобраны четыре структуры, отра- 
жающие структурный полиморфизм теломерной С-квадруплексной 
ДНК в виде четырех основных структурных топологий, а именно: 
ПРОВШ 1КЕТ, которая отражает топологию “ргореЙег”, 2) РОВШ 
2НУ9- топология “Вубма-1” 3) РОВ 2.$[ —“ы 9-2”, 4)РОВТО 143) 


-топология “ФазКеЁ” (Рис. 1). 











АССОТТАСОСТТАССОСТТАС ОВ ТАСОСТТАССОТТАССОТТАССОТТ АССОТТАСССТТАССОСТТАСОС  АААССОТТАСОСТТАССОТТАСОССАА 











Рис. 1. Структурные топологии теломерных С-квадруплексных ДНК. 
А) 1430 - “БазКеЁ”, длина последовательности 22н. 
В) 2.75 — “Шно-2>, (25 н); С) 1КЕ1 - “ргореПег” ( 22 н). 
О) 2НУ9 - “шЬча-1>, (26н). 
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Структуры теломерных С-ДНК были взяты из базы данных РОВ 
[9]. Обработка данных о структуре лигандов осуществлялась пос- 
редством программного пакета ОрепВабеу2.4.0. [10]. 

Для проведения виртуального скрининга и выявления сайтов 
связывания отобранных кандидатов с С-ДНК (молекулярный докинг) 
были использованны программные пакеты АчюБоскУша [11 и 
АчюросКТоо[$ [12]. Для визуалицаии данных использовались прог- 
раммные пакеты РУМОГ. [13] и УМО [14]. В качестве базы данных 
низкомолекулярных лигандов, стабилизирующих теломерную С-ДНК, 
была использована специализированная база данных МАГОВ [15]. 

Скрининг проводился в 2 этапа. На первом этапе суммарное ко- 
личество отобранных соединений составило 24341. На втором этапе 
был проведен повторный скрининг с учетом всех имеющихся конфор- 
маций лигандов 4 потенциальных кандидатов, продемонстривовавших 
наилучшие показатели на первом этапе, что в сумме представляет 386 


соединений. 


Результаты и обсуждение 


По итогам проведенного виртуального скрининга были выявле- 
ны соединения, связывающиеся с высокой аффинностью с рядом опи- 
санных теломерных С-квадруплексных структур (Табл. 1). 

Из выявленных соединений приоритетными являются те, кото- 
рые проявили высокую энергию связывания одновременно с различ- 
ными теломерными С-ДНК и допустимые показатели КМЗО. С уче- 
том вышесказанного, наибольший интерес представляет одно из про- 


изводных дикарбоксиамида* (Рис. 2). 
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Таблица 1. 


Результаты виртуального скрининга. Представлены топ 10 
потенциальных кандидатов для каждой структуры 
в порядке убывания энергии взаимодействия. 


Производное 
димиданафта- 
лена 





ВОО! 
(6408010) 


Производное 
триазола 


Производное 
нафталена 
(0408164) 


: Производноеди- 


Пиридил 
(04)0ВО6) 


: Производное ак- : 


ЗАОМ 






Производное 
дикарбоксиа- 
мида 


карбоксиами- 
да** 
(94080577) 


ВКАСО-19 
(04080285) 


- Производное ди- 


карбоксиамида 


: Производное ак- : 


ридина 
(0408280) 


РгореЙег : Нубла-1 
‘юроювзу (1КЕТ) : юро!ору (2НУ9) 
: Производное 
Теломестатин ‘ дикарбоксиами- | 
(04080377) да“ 
(94080578) 
обоев ОЛА 
т (048556) 
Е о 
(6408565) 04080166 
- Производное ди- ° Производное 


да* 


(44080556) 
карбоксиамида 


Производное Г- 
: пролина 
- (040804) 


Производное 
диэтанамина 
(04080804) 


Нубч9-2 


: В15-производное : 


карбазола 
(040Вр441) 


ВОО1 
(040810) 


ЗАОМ 


: - Производное Т.-- 
- дикарбоксиами- - 14 о 


пролина 
(040804) 


бАОМ 


ридина 


РРС-М-ТО 
(64080518) 


Производное 
дикарбоксиа- 
мида* 


РЬеп-0ОС3 
(064080648) 


: В15-производное : 


карбазола 


и 


Производное 


(64080574) 


РОС-М-ТО 
(64080518) 


: | Производное 
Е дикарбоксиамида: дикарбоксиами- 
да (904080560) 


- Производное ди- - 


миданафталена 
((40Вр27) 


Производное ди- 


карбоксиамида 
((40В0578)** 
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6АОМ а Теломестатин а 
(94080174) : (04080570) : (6408377) : (04680556 


Полные названия и данные о химических соединениях представлены в базе данных 
МАГОВ. 

Под соединениями указаны их идентификационные коды в базе данных. 
*-М№-В15(2-(рарегалт-1-У1)-асеёатл14о-6-руп9У1)-1,10-рвепапгойпе-2,9- 
Чсагбохапие; 
**-М,М№-В1$(2-(3-р1регахт-1-У])-ргорапап9о-6-рут4у1)-1,10-рБепапго[те-2,9- 


Фсагоохапие. 











[в 









































Рис. 2. Взаимодействие М,№'-В15 (2-(ррега7т-1-у-асеапи9д0-6-рутдУ)- 
1,10- рвепап го! пе-2,9-Ф1сагЬохапи4е с теломерными С-ДНК. 
А)143); энергия связывания (ЭС) = - 8.7 ккал/моль; 

В)25Л;; ЭС = -9.7 ккал/моль С) 2НУ9; ЭС = -9.5 ккал/моль; 
О)1КЕТ; ЭС = -9.3 ккал/моль. 

Показатели ВМЗ во всех случаях <= 0.2 нм. 
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Таким образом, в результате проделанных исследований был об- 
наружен ряд соединений, таких, как ЗАОМ, 6бАОМ и два производных 
дикарбоксиамида, взаимодействующих и стабилизирующих разные 
структурные топологии теломерной С-ДНК. Важно отметить, что в 
сравнении с лидирующим соединением ВКАСО-19, кандидаты проя- 
вили большую энергию связывания. При этом, ВКВАСО-19 свойствен- 
на высокая энергия связывания только по отношению к одной струк- 
туре теломерной С-ДНК, и в т ухо экспериментах соединение не 
продемонстрировало предсказанных свойств. Основная причина того, 
что ВВАСО-19 до настоящего времени считается лидирующим соеди- 
нением, связана с попыткой искусственно подавить влияние динами- 
ческого полиморфизма квадруплексных структур путем введения до- 
полнительных ко-растворителей типа ПЭГ, ацетонитрила, а также са- 
харов, белков, бивалентных катионов и т.д. с целью создания условий 
для формирования одной доминантной структуры. Такой подход в те- 
рапевтическом плане признан нецелесообразным [8]. Актуальность и 
важность настоящего и аналогичных исследований вне сомнений, т.к. 
позволяет предотвратить неоправданные и зачастую пагубные мани- 
пуляции, изменяющие микроокружение и биологию клеток и тканей. 

Работа выполнена при поддержке гранта МОН РФ Российско- 


Армянского университета (РАУ). 
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АННОТАЦИЯ 


В последние годы, следуя принятой тенденции использования раз- 
личных природных соединений, ведутся поиски ресурсов безвредных тера- 
певтических препаратов. Экспериментально показано, что цитрусовые фла- 
воноиды, включая полиметоксифлавоны, рассматриваются как потенциаль- 
ные химиотерапевтические агенты в борьбе с канцерогенезом и обладают 
широким спектром биологической активности, в том числе, ингибированием 
факторов риска метаболического синдрома. Целью данной работы являлась 
оценка возможного эффекта действия одного из самых распространенных 
полиметоксифлавонов — танжеретина на регуляторы апоптоза: гетеротетра- 
мер и гетеродимер белков Вс1-2 и ВНЗ домена ВАХ, а также на рецепторы 
смерти ОВК1,0ОЕ2, ОВ4, ОК5. Результаты показали, что танжеретин взаимо- 
действует с аминокислотами гетеротетрамера, участвующими в образовании 
гетеродимерных связей. Также показано взаимодействие с рецепторами 
смерти. Это дает возможность рассматривать полиметоксифлавон танжере- 
тин в качестве перспективного антиканцерогенного агента. 

Ключевые слова: танжеретин, регулятор, рецептор, апоптоз, докинг, 


аффиность. 
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Введение 


В последние годы, следуя принятой тенденции использования 
различных природных соединений, ведутся поиски ресурсов безвред- 
ных терапевтических препаратов. Потенциальное протекторное дейст- 
вие фруктов и овощей на здоровье человека связано с наличием в них 
вторичных метаболитов, биоактивных функциональных соединений, 
таких, как: флавоноиды, каротиноиды, аскорбиновая кислота и лимо- 
ноиды, обладающие антиоксидантной, противовоспалительной и про- 
тивораковой и др. активностями [1, 2]. Экспериментально показано, 
что флавоноиды цитрусовых могут действовать как профилактические 
и химиотерапевтические агенты самостоятельно, а также в качестве 
адъювантов с другими лекарствами [2]. Полиметоксифлавоны это 
группа метоксилированных фенольных соединений, встречающихся в 
основном в цитрусовых, которые состоят из бензо-гамма-пирона с 
карбонильной группой в положении С-3 и метоксигрупп в разных по- 
ложениях. Биологическое значение полиметоксифлавонов, по сравне- 
нию с другими флавоноидами, выше благодаря липофильным свойст- 
вам большинства находящихся в них метоксильных групп [3]. Флаво- 
ноиды цитрусовых рассматриваются как потенциальные химиотера- 
певтические агенты в борьбе с канцерогенезом и обладают широким 
спектром биологической активности, в том числе, ингибированием 
факторов риска метаболического синдрома [4—6]. 

Целью данной работы являлась оценка возможного эффекта 
действия одного из самых распространенных полиметоксифлавонов- 
танжеретина на регуляторы апоптоза: гетеротетрамер и гетеродимер- 
белков В<1-2 и ВНЗ домена ВАХ, а также на рецепторы смерти ОВ1, 
ОБ2, ОВ4, ОВ5. 
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Материалы и методы 


Для проведения докинга использовалась программа Ашюдоск 
Уша[7], используящая ламарковский генетический алгоритм 
(Гатайчап Сепейс А|еогит, ГСА). При создании модели молекулы 
танжеретина были получены с использованием разных программ по 
биомоделированию, а также банки данных КСЗВ Ргоеш Ба ВапК [8], 
РибСВет [9]. Структура рецептора была взята из банка данных ВСВ 
РОВ (Ш 2ХАОД). 

В качестве базовых параметров для гетеротетрамерной структу- 
ры были выбраны виртуальный бокс с размерами 76Х50Х56 и конти- 


нум исчислений 128. 


Результаты и обсуждение 


В работе проведен молекулярный докинг танжеретина с гетеро- 
тетрамером гетеродимера белка Вс!-2 и ВНЗ домена белка ВАХ, 
состоящим из четырех цепей А, В, С, О. Цепи А и В идентичны и при- 
надлежат анти-апоптотическому белку Вс|-2, С и О тоже идентичны и 
принадлежат ВНЗ домену (52-82а.о.) про-апоптотического белка ВАХ. 
Результаты докинга показали, что танжеретин взаимодействует с ами- 
нокислотными остатками [ГЕО59, СГМ118, НГ$120, АВО139 гетеро- 
тетрамера (Рис. 1). Из них ГЕО59 локализован в ВНЗ (а2) домене бел- 
ка ВАХ, (0[М118 в ВНЗ домене, а Н!$120, АКС139 в ВН! домене бел- 
ка Вс|-2. 

Проведенный нами докинг позволил получить модель взаимо- 
действия танжеретина с гетеротетрамером белков В‹с|-2 и ВНЗ домена 


ВАХ со средней энергией связывания -6.0- 0.216Ксало| 
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6) 


Рис. 1. (а) Докинг танжеретина в гетеротетрамере 2(ВНЗ ВАХ -Ве-2) 
А-В цепи Ве]1-2 и (С-В цепи ВАХ). 
(6) Аминокислотные остатки гетеротетрамера, вовлеченные 


во взаимодействие. 
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При анализе связывания танжеретина с гетеротетрамером В‹|-2 и 
ВАХ докинг имеет ряд недостатков, результаты менее надежны. Для 
уточнения правильности места связывания был проведен докинг с ге- 
теродимером, состоящим из А цепи белка Вс]1-2, и С цепи белка ВАХ. 
Энергия связывания, при этом, равнялась -5.6-0,216 Ксато]. Это сви- 
детельствует о том, что взаимодействие танжеретина с гетеротетраме- 
ром сильнее, чем с гетеродимером Вс|-2 и ВАХ (Рис. 2, а). 

При этом, танжеретин электростатически взаймодействует с ами- 
нокислотными остатками УАГ, 133, [/У558, ОГУ128, АВО129. Из них 
Г.У$58 локализован в ВНЗ домене белка ВАХ(С-цепь), а УАГ 133, 
01128, АВС129 в А-цепи белка Вс|-2 (Рис.2, 6). 


в 





Рис. 2. Взаимодействие танжеретина с гетеродимером 
ВНЗ ВАХ -Ве]-2. А цепь Вс]-2-тёмно-серый цвет и С цепь ВАХ-серый (а). 


Аминокислоты гетеродимера, вовлеченные во взаимодействие (6). 


Сравнение аминокислотных остатков гетеротетрамера и гетеро- 
димера, вовлеченных во взаимодействие с танжеретином, показывает 


на участие расположенных по соседству аминокислот (Табл.1). Танже- 
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ретин при докингес гетеротетрамером взаимодействует с АКС139 и 
ОТО 118 цепи А, которые образуют водородные связи с аминокисло- 
тами цепи С, тогда как при докинге в гетеродимере взаимодействует с 
ОГО 128 и АКС 129, которые участвуют в образовании водородных 
связей гомодимера А-Б, [1$ 58 цепи С, также вовлеченным в образо- 
вание гомодимера С-Д и лишь УАГ. 133 вовлечен в образование гете- 
родимера. Таким образом, танжеретин взаимодействует с аминокисло- 
тами одного и того же сайта, но электростатические взаимодействия с 
гетеротетрамером сильнее и комплекс стабильнее, а энергия взаимо- 
дейтвия с выше, чем с гетеродимером (-6.0- 0.216КсаМто] Ксалто! >- 
5.6+0.216 Ксато]). 

Таблица 1. 

Аминокислотные остатки гетеротетрамера и гетеродимера, 


вовлеченные во взаимодействие с танжеретином. 























Мишени Аминокислоты Аффинность 
(ккаль/моль) 
2хХВАХ( ВНЗо2) | Лейцин 59(а2) 
гетеротетра- Ре 
мер Глутамин 118(@3- а 4) 6.0+ 0.216 
2хВСГ-2 Аргинин 139 94- 9 5 
Гистидин 12093- 94 
ВАХ( ВНЗо2) Лизин 58 (42) 
гетеродимер В 
Глицин 128 04 -5.6=0.216 
ВСГ-2 Аргинин 129 94 
Валин 133 4 




















Взаимодействие с гетеротетрамером затрагивает аминокислот- 
ные остатки белка Вс|-2 а2 спирали домена ВНЗ и бета повороты меж- 


ду 03- а4 и 94 -а 5 спиралей домена ВНТ, с одной стороны, и ГЕЦ 59 
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а2 спирали домена ВНЗ белка Вах, вовлеченных в образование гетеро- 
тетрамера [10]. 

Это может дестабилизировать гетеротетрамер и способствовать 
удалению Вах, @2 спирали которого сосредоточены в области 59-101 
а.о. функционального ВНЗ домена, который достаточен как для гомо- 
димеризации, гетеродимеризации с белками семейства Вс|-2 [11], так 
и для стимуляции апоптоза в клетках млекопитающих. 

При взаимодействии с гетеродимером танжеретин, в основном, 
затрагивает вовлеченные в гомодимеризацию аминокислоты, где, в ос- 
новном, задействованы аминокислотные остатки спирали 94. 

При апоптозе М№-конец Вах-а подвергается конформационным из- 
менениям, которые затрагивают а@- 5 и 0- 6, которые могут быть су- 
щественны для выворачивания белка наизнанку и включены в актива- 
цию белка, сопровождающуюся ассоциацией и образованием пор на 
митохондриях. Основываясь на сильных молекулярных контактах 
между спиралями 0- 1, 2, 5, 6, и 9, изменения в конформации любого 
из этих регионов теоретически должны триггерировать изменения в 
других регионах. В частности, поскольку 092, содержащий ВНЗ домен, 
примыкает к обоим терминальным 01 и 99 спиралям, то смещение лю- 
бого из них может вызывать экспозицию ВНЗ. 

Апоптоз представляет собой высокозапрограммированную кле- 
точную смерть для ремоделирования ткани, гомеостаза и развития 
многоклеточных организмов. Одним из целевых сайтов для лечения 
канцерогенеза является активация внешнего апоптотического пути, 
т.к. при этом клетки неопластической ткани становятся нечувстви- 
тельными к апоптозу в зависимости от гистопатологии, но, чаще все- 
го, у них развивается резистентность к СО95 (АРО-1, Еа$) апоптозу. 
Внешние пути апоптоза связаны с кластеризацией рецепторов смерти 
на поверхности клетки после активации взаимодействием с тримерами 


лигандов. Все рецепторы смерти принадлежат суперсемейству рецеп- 
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торов фактора некроза опухоли (ТМЕК): ТМЕ -В1, СО95 (АРО-1, Еа$), 
ТВАМР (АРО-3 ОК-3), рецептор 1, связанный с ТМЕ апоптоз-индуци- 
рующим лигандом (ТВАП)) - (ТВАП.-В1 ОВ4) и ТВАП--В2 (АРО-2, 
ОК5). Для достижения лучшего терапевтического действия лидерной 
молекулы, а также уменьшения побочных эффектов молекул в других 
регионах предпочтительной мишенью являются рецепторы смерти. 
Поэтому нами был проведен докинг-анализ танжеретина с рецептора- 
ми смерти ОВ1, ОК2, 284, ОК5. В Табл 2. обобщены результаты до- 
кинг-анализа танжеретина с аминокислотными остатками ОК 1[1ЕЗУ]; 
ОК 2- [3МХ7]; ОБ4- [10206] ОК 5- [104У]. 


Таблица 2. 
Докинг танжеретина с аминокислотными остатками ОВ 
Т1ЕЗУ|; ОВ 2- [3МХ7|; ОВ4- 106] ОВ 5- [104У|. 












































Лиганд | Мишень Аминокисло | Онаьох А Средняя энергия 

ты взаимодействия 
Танже- | МТДТВАРР |01[0151 76*44*54  |-5,815= 0.318Кса/тпо| 
ретин и ТКАЕ домен 

ТВАЕ2 

Танже- | ЕАМ-МТО ЗЕВ85, 110*76*96 |-4,72-=0,179Ксато1 
ретин ГЕЧ12, 

ГУ5$84 
Танже- | АРО2Г/ТКАП, | АКО121, 106*80*101 |-5,38=0,150Ксато1 
ретин АГА123, 

РНЕ278 
Танже- | ТКАП, -ОВ5 ТУБ189, 118*90*118 |-5,625-0, 183КсаИло1 
ретин АГА272, 

©1271 

Как видно из табл. 2, танжеретин электростатически взаймо- 





действует с аминокислотным остатком ©1151 ОВ1; $ЕК85, [ЕО12, 
ГУ$84 — ЕАШМ-МТО; АКО121, АГА123, РНЕ278 АРО2Г/ТВАЦП,, а 
также ТУК 189, АГ.А272, 91271 ТВАП, -ОВ5. 
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Таким образом, полученные результаты свидетельствуют о том, 
что танжеретин взаимодействует с аминокислотами гетеротетрамера 
Вс1-2 -ВНЗ ВАХ, участвующими в образовании гетеродимерных свя- 
зей, тем самым, способствуя освобождению про-апоптотического бел- 
ка. Также показано взаимодействие с рецепторами смерти ОК1, ОВ2, 
ОБА, ОВ5, которые участвуют в регуляции апоптоза. Это дает возмож- 
ность рассматривать полиметоксифлавон танжеретин в качестве перс- 


пективного антиканцерогенного агента. 
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п гесепЕ уеагз, гоПо\уте Ше тепа оЁ Фе и$е оЁ уапоц$ пабага] сотпрочпа$, 


фе зеагсВ г Вагп[ез$ Фегареиис 4гиез 1$ Беше сопаис{е4. И Ваз Бееп зпо\п 
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шуо[уе4 ш бе югтайоп оЁ Веегодитегс Боп4$. ПиегасНоп \ИВ деай гесерюг$ 
15 а15о зпо\т. ТЬ1$ таКез # роз$1Ые {ю феа{ ро]утеохуНауопе {апоегейп аз а 
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АННОТАЦИЯ 


Артемизинины проявляют широкий спектр биологической активности, 


такой как антиоксидантная, антипротозойная, антиканцерогенная, противо- 


грибковая, противовирусная, противоспалительная и т.д. Однако до настоя- 


щего времени механизмы действия и молекулярные мишени артемизинино- 


вых недостаточно изучены. Предварительный виртуальный скрининг биоми- 


шеней артемизинина выявил сильное взаимодействие последнего с лиганд 


связывающим доменом глюкокортикоидного рецептора (ВО СВ) человека, 


с которым были проведены докинг, анализ по главным компонентам (РСА), 


кластерный анализ и оценка кластеризации. Результаты показали, что арте- 


мизинин связывается с (ВО СВ в трех сайтах, среди которых наибольшее 


значение показали первый и третий (65% и 20%, соответственно) и затраги- 


вает чрезвычайно важные сайты [ВО СВ. 


Ключевые слова: артемизинин, ГВО глюкокортикоидного рецептора, 


докинг, РСА, кластерный анализ. 
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Введение 

Аптептияа аппиа Г., широко известная как сладкая полынь или 
ОтеВао, крупный ядовитый однолетний кустарник листья которого 
выделяют эфирные масла, используемые в народной и современной 
медицине, а также косметической и фармацевтической промышлен- 
ности. Некоторые виды полыни, включая А. аппиа, традиционно ис- 
пользовались в Китае при лечении 52 болезней, болях, воспалениях и 
лихорадочных состояниях. Основными действующими компонентами 
полыни однолетней являются вторичные метаболиты артемизинин и 
его производные [1]. Артемизинины принадлежат семейству сескви- 
терпеновых триоксановых лактонов, проявляют антиоксидантный, 
противовоспалительный, антиканцерогенный, иммуномодулирующий, 
противомикробный, противогрибковый, антигельминтный, противови- 
русный и другие активности. Однако до настоящего времени механиз- 
мы действия и молекулярные мишени артемизининовых недостаточно 
изучены [2]. Предварительно проведенный нами скрининг биомише- 
ней артемизинина выявил сильное взаимодействие последнего с ли- 
ганд-связывающим доменом глюкокортикоидного рецептора (СВО 
ОК) человека (РОВ:49ОС). 

ОВ является конститутивно экспрессируемым транскрипцион- 
ным фактором, который контролирует множество различных генных 
сетей, каждая из которых уникально детерминируется конкретным 
клеточным и физиологическим контекстом [3]. 

Целью представленной работы было т 5Шсо исследование взаи- 
модействия артемизинина с ГВО СК методами молекулярного моде- 


лирования. 


Материалы и методы 
При проведении докинга использовались программные пакеты 


[4]. Анализ по главным компонентам (РСА) параметров конформации 
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результатов докинга проводили по методу [5], кластерный анализ — ал- 
горитмом К-средних [6]. Качество кластеризации оценивалось индек- 
сами Дэвиса-Болдина (ОВ), Данна (Пипп-шдех), оценки силуэта 
(ЗПРочейне шдех) и Калински-Харабаза (СН) [7]. Для выявления сайтов 
связывания артемизинина на мишени была использована программа 
ГОРГОТ [8], для визуализации данных — РУМОГ [9]. 


Результаты и обсуждение 

Результаты анализа докинга и кластеризации показали, что арте- 
мизинин связывается с [ВО СК в трех сайтах, из которых большая 
часть принадлежит первому кластеру (65%) (рис. 1, а). При этом энер- 
гия связывания артемизинина с ним наивысшая. Оценка кластеризации 
выявила, что значение индекса Данна низкое, а оценки индексов Дэви- 
са-Болдина, силуэта и Калински-Харабаза (рис. 1, 6) — довольно высо- 


кие, что свидетельствует об оптимальном качестве кластеризации. 
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Рис. 1. а) РСА с кластеризацией на основе данных докинга 
с использованием алгоритма К-средних. 


6) Определение оптимального количества кластеров. 
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Рецепторы всех стероидных гормонов имеют общую доменную 
архитектуру, состоящую из иммуногенного М-концевого домена 
(МТО), ДНК-связывающего домена с цинковыми пальцами (ОВО), 
шарнирной области и С-концевого лиганд-связывающего домена 
(ГВО), который состоит из 12 альфа-спиралей (Н) и 4 В-складок, обра- 
зующих гидрофобный карман для связывания лиганда [3]. (ВО также 
содержит спираль АЕ-2, проявляющую функцию лиганд-зависимой 
активации, регулируемой, например, связыванием с гормоном [10]. 
Кластерный анализ молекулярного докинга артемизинина с ГВО СК 
(таблица 1), а также анализ при помощи программы ГОРГОТ пока- 
зали, что в первом сайте А, В, С и В кольца артемизинина образуют 
гидрофобные связи с а — спиралями НЗ, Н7, Н10/Н11 и Ва-листом [ВО 
ОВ. 


Таблица 1. 
Энергии связывания (Ккал/моль) кластеров, 


содержание каждого в %, аминокислоты [ВО СВ 


и вовлеченные во взаимодействие кольца артемизинина. 





Клас- | Энергия свя- | Содер- Аминокислотные Кольца ар- 
тер зывания жание в остатки [ВО СК. темизинина 
(Ккал/моль) % 





1 -7,069+0,479 65 Ме(560 (НЗ), Геи563 А, В, С, О 
(НЗ), Азп564 (НЗ), 
РНеб23 (ВА), 911642 
(НТ), Меб46 (НТ), 
[.еи732 (Н1О/НИ1), 
Туг735 (Н10/Н!\), 
Пе747, Рве749 
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2 -6,567+0,287 5 Гуз695 (Н9), ©1уб98 А, С, Б 
(Н9), Пе701 (Н9), 
\а1702 (Н9), Азп711 
(Н1О), РБе7 15 (Н10) 
3 -6,2+0,071 20 ©1542, \Уа1543, Тгрб10 | А, В, О 
(Н5), Аг2611 (Н5), 
Атс614 (Н5), Тугб60 
(Н8) 
































*данные для 2000 конформаций. 


Известно, что спираль НЗ и Вл-лист содержат сайты сильной диме- 
ризации [ВР СК, а аминокислоты спирали Н5 участвуют в стабилизации 
димера [10]. В участке НЗ-Н7 находится минимальный сегмент (550-653 
а.о.), необходимый для высокоаффинного связывания с шапероном 
Нзр90 [11]. В спиралях НЗ, Н5 и Н!2 расположены сайты связывания с 
коактиватором М№соА-2 [3]. Следовательно, взаимодействие артемизини- 
на с СК должно отражаться как на димеризации и стабилизации Г.ВО, так 
и на связывании с Нзр90 и коактиватором №соА-2. Артемизинин взаимо- 
действует с аминокислотными остатками 1747 и Е749 петли, предшест- 
вующей спирали АЕ-2. Известно, что это сопровождается стабилизацией 
спирали АЕ-2 в активном состоянии и может представлять молекуляр- 
ную основу для лиганд-зависимой активации ОК [10]. 

Во второй сайт гидрофобных взаимодействий вовлечены А, С, О 
кольца артемизинина и аминикислоты спиралей Н9 и Н1О (таблица 1). 
В третьем сайте Тгрб10, Аго611, Аго614 спирали Н5 и Тугб60 спирали 
Н8 образуют гидрофобные связи с кольцами А, В, О (таблица 1). Учи- 
тывая роль спирали Н5, это приведет к стабилизации димеризации 
ЕВРО СК. 

Эндопероксидный компонент в структуре артемизинина (С коль- 
цо) играет важную роль в его биологической активности [12]. Взаимо- 


действие с первым и вторым сайтами ГВО СВ также опосредуется С 
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кольцом, возможно как из-за его реакционноспособности, так и спо- 
собности к димеризации. 

Артемизинин приводит к экспрессии конститутивного андроста- 
нового рецептора (САК) и белков семейства цитохромов Р450 (СУР5), 
которые играют решающую роль в регуляции метаболизма многих ле- 
карств, ксенобиотиков, включая артемизинины, гомеостаза энергии и 
неопластического перерождения модулированием транскрипции мно- 
гочисленных генов-мишеней [13]. СУРЗА4/5, СУР2С9, СУР2С19 и 
СУР2О6 являются клинически значимыми для метаболизма лекарств 
[14]. С другой стороны, активация СВ сопровождается как индукцией 
САК, так и СУРЗА. Обобщая полученные нами результаты и со- 
поставляя с имеющимися данными по механизму действия артемизи- 
нина, нами предложена гипотетическая схема взаимодействия СВ-ар- 
темизинин с энхансерными элементами промотора СУРЗА5 с экспрес- 


сией гена (Рис. 2). 





Рис. 2. Гипотетическая схема 
| действия артемизинина на 
экспрессию гена СУРЗА$, 
опосредованного связыванием 
.ВО СВ с энхансером 
промотора (-891/-1109, два полу- 


сайта глюкокортикоид- 





зависимых элементов (СВЕ)). 


-891/-1109 (два полусайта) 
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Таким образом, нами впервые методами молекулярного модели- 
рования показано, что артемизинин связывается с ГВО СК в трех сай- 
тах, среди которых наибольшее значение показали первый и третий 
(65% и 20%, соответственно) и затрагивает чрезвычайно важные сай- 
ты ГВО СК. На основании полученных результатов можно прийти к 
выводу, что действие артемизинина на множественные сигнальные пу- 
ти клетки может быть опосредовано глюкокортикоидным рецептором, 
а артемизинин можно рассматривать в качестве нового лиганда для 


глюкокортикоидного рецептора. 
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АННОТАЦИЯ 


Вопросы систематики мышей рода Аро4етиу зепзи 1афо остаются до 
настоящего времени одними из самых сложных в современной таксономии 
грызунов. Несмотря на многочисленные работы, таксономический статус, в 
особенности, биогеография рода до сих пор не уточнена в Закавказье, в част- 
ности, на территории Армении. В данной статье приведен анализ нуклеотид- 
ных последовательностей с построением медианной сети гаплотипов степной 
мыши Аройетих уййегЬу для уточнения таксономического статуса вида на 
территории Армении. 

Ключевые слова: таксономия, лесные мыши, Ародетиз, биоразнооб- 


разие, Армения. 


Введение 


Группа лесных мышей представляет собой многочисленную и 
полиморфную группу мышевидных грызунов, населяющих различные 
ландшафтные пояса республики Армения. Лесные мыши рода 


Аро4етиб зепзи ао широко распространены на всей территории Ар- 


В.Б. Азарян 145 





мении. Информация о таксономии и распространении видов рода 
Ародетиз на территории Армении является весьма скудной. По пос- 
ледним данным исследователи фауны Армении выделяют 5 видов лес- 
ных мышей, а именно — Аро4етиу игепя5, Аро4етиу илйегБуй, 
Аро4детиб ропйсиз, Аро4етиу ПаялсоШ5, Аро4етиз 5угайсиз. 

И.Я. Павлинов с соавт. [1-2] придают статус рода двум подродам 
группы лесных и полевых мышей — Ародетиз и Зууаети$, а 
С.В. Межжерин [3-4] — трем (Аро4етиз, Зулаетиз, Аботуз). Западные 
систематики, в частности, Миззег, Сайеюп, Егуша е{ а|., придержи- 
ваются прежних взглядов о единстве рода Аройетиб зепзи [а1ю, выде- 
ляя в его составе лишь подроды или группы видов [5, 6-7]. Первона- 
чально все виды лесных и полевых мышей причислялись к единому 
роду Аро4етиз, а в наши дни на этот счет существуют серьезные воз- 
ражения. Мы придерживаемся точки зрения Егутиа её а[., используя 
название политипического рода Аро4етиз зепзи 1а0 как для группы 
западно-палеарктических лесных мышей булаетиз (с подродами 
5у’аетиб и Кагуюту5), так и для восточно-палеарктических лесных и 
полевых мышей Аройетих зепза 5И/ (с подродами Аро4етиз и 
Ароту5). Ситуация осложняется неравномерной изученностью видов 
данной группы, а также продолжающимся первоописанием новых 
форм [8—10, 11]. Ученые не отрицают также мнение о возможности 
нового гибридогенного видообразования [12—15], как и продолжения 
формообразовательных процессов в популяциях лесных мышей и по 
сей день [10, 11, 16]. Несмотря на то, что систематика лесных и поле- 
вых мышей еще далека от завершения, тем не менее названия по всему 
тексту приведены в соответствии с современной западной зоологичес- 
кой номенклатурой: Аро4етиз зепзи 1афо. 

В данной работе представлены новые генетические данные по 
лесным мышам рода Аро4ети5 в Армении. Полученные путем секве- 


нирования митохондриального гена первой субъединицы цитохромок- 
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сидазы (СОТ) нуклеотидные последовательности были использованы 
для конструирования сети гаплотипов вида Аройетиз жийегБуг. Целью 
данной работы являлось: 1) уточненная видовая идентификация соб- 
ранных образцов вида молекулярно-генетическим методом; 2) анализ 
нуклеотидных последовательностей на популяционном уровне; 3) об- 


суждение полученных результатов. 


Материалы и методы 


Было проанализировано в общем количестве 77 образцов мышей, 
принадлежащих к виду степной мыши Аройетих у’ИйегЬуг. Географи- 
ческая принадлежность образцов представлена следующими локалите- 
тами (Таблица 1): 1. Агавнадзор (40°35'03"М 44°41'29"Е; количество об- 
разцов п = 10), Котайская область, Центральная Армения, 1900 м над 
уровнем моря. 2. Айгут (40°41'03"М 45°10'23"Е; п = 25), Гегаркуникская 
область, Северо-Восточная Армения, 1400 м над уровнем моря. 3. Ере- 
ван (40°11”М 44°3 ТЕ; п = 3), Центральная Армения. 4. Хосровский на- 
циональный заповедник (39°56'42.1"М 44°51'35.5"; п = 5), Араратская 
область, Центральная-южная Армения, 5. Шванидзор (38°54’10”М 
46°14'40”"Е; п = 24), Сюникская область, Южная Армения, 657 м над 
уровнем моря, 6. Ванадзор ( 40°48'16"М 44229'38"Е; п = 10), Лорийская 
область, Северо-Восточная Армения, 1344 м над уровнем моря. 


Таблица 1. 


Список локалитетов генетически определенных образцов 











Ародети$ уупегру!. 
А. уйдетБу 
Количество Количество 
экземпляров (М) гаплотипов (№) 
Котайкская область 10 5 
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А. уфетБу 
Количество Количество 
экземпляров (М) гаплотипов (№) 

Гегаркуникская область 25 12 
Ереван 

Араратская область о 5 
Сюникская область 24 13 
Лорийская область 10 5 

















Общая геномная ДНК была изолирована из образцов хвоста или 
ногтя, зафиксированных в абсолютном этаноле при — 20°С. Изоляция 
ДНК была проведена с использованием ОМАеазу В]оо4 апа Т15зе КИ 
(О1асеп), следуя протоколу производителя. Амплификация ДНК была 
произведена ПЦР реакцией, используя специальные праймеры для 
мтДНК. Все ПЦР-продукты были очищены осаждением в этаноле. По- 
лученные нуклеотидные последовательности (секвенсы) были обрабо- 
таны и выравнены вручную, используя программу В1оЕай. Выравнен- 
ная нуклеотидная последовательность, имеющая длину 728 п.о., была 
подготовлена для реконструкции сети гаплотидов. Был использован 
метод теФап-]отше, доступный в программе МЕТ\У!ОВК, уегзоп 
4.6.1.2 для конструкции сети гаплотипов анализируемого вида 
Аро4етих ийетБуг. Медианная сеть является графическим представ- 
лением последовательности мутационных шагов, в результате кото- 
рых общая предковая последовательность нуклеотидов превратилась в 
совокупность родственных последовательностей ныне живущих по- 


ТОМКОВ. 


Результаты и обсуждение 


Для вида Ародетиз уййе’Бут была сконструирована выровненная 


нуклеотидная последовательность, состоящая из 728 пар оснований, 
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из которых 31 были вариабельными и 55 экономно-информативными. 
Данная выравненная последовательность была использована для 
конструирования сети гаплотипов (Рис. 1) и расчёта характеристик 
гаплотипного разнообразия. В исследованных выборках было обнару- 
жено 43 гаплотипа для 77 особей. В популяциях лесных мышей отме- 
чено характерное распределение гаплотипов, т.е. небольшое их число 


имело высокую частоту, другие варианты встречались редко. 


вы» — Котайская область 
Араратская область 


Сюникская область 
= Баш :г^ 


Лорийская область 
Гегаркуникская область 45 © 


$355. 
А О у © - 385 
\ 793АР® 
| о ик. 
1 ®] г =) 
влзир ® - ао 
а ЗА те 5^@ 
мал © кА 
еж 7323АР Эм зал @} | 7ОЗАР 
тж 
ПАРФ зи 

о зи 


Рис. 1. Медианная сеть гаплотипов вида /4. жийегруг, полученная 
для последовательностей митохондриального гена первой 


субъединицы цитохромоксидазы (СОТ). 


Степная мышь — самый сухостойкий и, очевидно, самый тепло- 
любивый вид рода лесных мышей. Вид Аро4етиз уийетЬут (ТБотаз, 
1902) был описан в Израеле [11, 17] как А. йегтопеп515 [11], однако, 
позже оказалось, что он широко распространен по степным зонам 
Среднего Востока (сте ВВодо$ 151апа; [18]). Позже Филиппуччи с 
соавторами (1996) [19] предложили, что А. йегтопеп15 скорее всего 


может быть конспецифичным с видом ранее заявленным как А. 
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лает (Ме77Вегт & Гасогодпучик, 1989) или А. /и[яресиу (Озпеу, 
1924) из Туркменистана, Закавказья, Кавказа и соседних степей к Кры- 
му [20]. Наличие А. Регтопеп515 в Армении было заявлено Сузуки с 
соавт. (2008) [21]. В настоящее время западные исследователи счи- 
тают А. йегтопеп515 малым синонимом к Аройетих илйетБуг [22, 23, 
24, 25]. Ареал степной мыши (Межжерин, 1997) охватывает Предкав- 
казье (вплоть до Прикаспийских степей и Калмыкии), Кавказ, Закав- 
казье› Малую и Переднюю Азию (Израиль), Северный Иран, пред- 
горья Копетдага. Внешние характеристики образцов, входящих в гап- 
логруппы, частично совпадают. Общая длина тела 94,3 мм; длина уха 
14,5 мм; длина задних лап 21,2 мм; длина хвоста 100 мм. У вида наб- 
людается бледная и более желтая дорзальная поверхность. Внешние 
размеры и окраска наших генетически определенных образцов А. 
уийет’Бут соответствуют вариации, описанной Ктубек & УоргайК 
(2009) [24]. 

Медианная сеть гаплотипов, выявленная с помощью анализа 
нуклеотидных последовательностей, даёт нам повод предположить не- 
давнее демографическое расширение популяций. Наличие некоторых 
достаточно своеобразных гаплотипов у исследованного вида требует 
дальнейшего детального изучения. География и климатическая исто- 
рия Армении позволяют предположить, что эта картина может быть 
результатом последовательных успешных волн колонизации. 

В работе была проведена не только морфологическая (повсе- 
местно используемая методика определения видов на территории рес- 
публики), но и генетическая идентификация вида. В общем и целом 
можно с большой вероятностью утверждать наличие и самостоятель- 
ность представленного вида на территории Армении, принимая во 
внимание представленный обзор литературы по данному вопросу и ос- 
новываясь на собственных результатах молекулярно-генетических 


исследований собранного материала. 
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АННОТАЦИЯ 


Формирование синантропизма — процесс сложный и многоэтапный. 
Млекопитающих, которые образовывают синантропные популяции, с каж- 
дым годом все больше. Однако, у разных отрядов этот процесс протекает по- 
разному. Это связано с образом жизни, питанием животных, а также с раз- 
личными антропогенными факторами. В данной статье опубликованы ре- 
зультаты наших исследований по трем отрядам — Насекомоядным 
([тзесйуота), Рукокрылым (Сйторега) и Грызунам (Ко4депна), приведен спи- 
сок основных синантропных представителей каждого отряда. 

Ключевые слова: синантропизм, антропогенные факторы, хищники, 


грызуны. 


Процесс формирования синантропизма у представителей разных 
отрядов млекопитающих изучен недостаточно. Именно поэтому 
целью данной работы является изучение видового состава синантроп- 
ных млекопитающих и их экологии в Армении и НКР. В нашей работе 


представлены материалы по наблюдению за животными, которые 
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были собраны с помощью фотоловушек. Приватизация земель и их 
бесконтрольная обработка, вырубка лесных массивов, выпас скота, а 
также преследование и обстрел хищных видов зверей во многих райо- 
нах Армении и НКР привели к тому, что некоторые виды животных, 
спасаясь от губительной деятельности человека, были вынуждены 
покинуть природные ландшафты, переселиться ближе к человеку, 
нередко, даже под одну крышу с ним. 

Люди с целью подготовки приобретенных под пашни участков 
земель, выжигают на них кустарниковую и травяную растительность, 
а затем, после сбора урожая, поджигают и остатки зерновых, якобы 
тем самым удобряя почву и уничтожая всевозможных вредителей. Эти 
процессы способствуют переходу с этих мест под укрытие всевозмож- 
ных ближайших строений человека не только грызунов и насекомояд- 
ных, но и хищников, рептилий и многих видов других членистоногих. 

Вырубка крупных дуплистых деревьев вынуждает рукокрылых 
(Сшторлега), многие представители которых обычно живут, размно- 
жаются, а нередко и зимуют под корой или в дуплах, искать убежище в 
чердаках домов, в подвалах, в сараях, овчарнях или за наличниками 
окон и дверей. 

Как это не парадоксально, всё совершенствующиеся методы и 
средства борьбы со всевозможной живностью, стремящихся «под кры- 
ло» человека, до сих пор не дали ожидаемого результата. 

Так, если всего несколько десятилетий тому назад человек до- 
вольно успешно, самыми примитивными средствами в своём доме и 
ближайшем окружении боролся всего с несколькими видами грызунов 
(крысы, мыши, хомяки), а насекомых изводил соляркой и керосином 
или щелочным раствором и селитрой, то сейчас не спасают даже спе- 
циальные пластиковые трубки от насекомоядных, живоловки, всевоз- 
можные давилки, капканы, а разбрызгиваемые аэрозольные средства 
нередко опасны для самого человека и его домашних животных (птиц, 


кроликов, голубей, кошек и собак). 
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Помимо вышеописанной хозяйственной деятельности, люди 
строят многоэтажные, панельные, многощелевые дома, в которых 
функционируют мусоропроводы — накопители отбросов. Все эти 
«удобства», без учета желания человека, являются прекрасным и удоб- 
ным местом для жилья, питания и размножения для большого числа 
видов летучих мышей, крыс, мышей, землероек и насекомых. В ре- 
зультате, имеем неожиданный и, конечно, нежелательный результат — 
рост синантропных видов. Причем, если некоторые из них в природ- 
ных условиях размножались лишь в определённое (теплое) время года, 
то перейдя в «тепличные» условия к человеку эти же виды не теряют 
активности и плодятся круглый год. 

Наши исследования в Армении и НКР в этой области, начатые в 
1991 году (Авагян, 2010; Попов, 2003; Явруян, 1991; Явруян, Арутю- 
нян, 1997; Явруян, Арутюнян, Казарян, Авагян 2006; Явруян, Айрапе- 
тян, 2003) и продолжающиеся по сей день, позволяют представить 
список современных основных видов млекопитающих-синантропов. 


Данный список представлен в Таблице 1. 
Таблица 1. 


Основные синантропные представители по трем отрядам. 








лая бурозубка (5отех 





ттиши$) 








Отряд Насекомояд- | Отряд Рукокрылые Отряд Грызуны 
ные (175есйуога) (Сштгориега) (Ко4епийа) 

Бурозубка Волнухи- | Большой подковонос | Соня-полчок 

на (5огех (Н.) (Кйтоюрйи5 (Муохиб (°.) 215) 

Готисйт)) или ма- |/етгитедитит) 
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Отряд Насекомояд- | Отряд Рукокрылые Отряд Грызуны 
ные (]пизесйуога) (Сшториега) (Ко4епиа) 
Кавказская бурозуб- | Малый подковонос Лесная соня 
ка (5огех (5.) (Кйтоюрйи5 (Огуотуз ппеаша) 
5аитт) - найдены 2 | ррозщето5) 


особи в подвале жи- 
лого дома в Айгуте. 





Остроухая ночница 
(Муой5 Буй) 


Серая крыса (КаНи5 
погуее1си5) 





Малая белозубка 
(Стосщига 
5иауео[еп5) 


Водяная ночница 


(Муой5 ааиБетопй) 


Чёрная крыса 
(Вагии$ гаци$) 





Белобрюхая бело- 
зубка (Сгосчига 
[еисодоп) 


Ночница Наттерера 
(Муой5 паЦеге!") 


Домовая мышь (Ми5 
тизси[и$) 





Длиннохвостая бе- 
лозубка (Сгосига 
гие4епяае4!") 


Трехцветная ночница 
(Муой5 етагета!и5) 


Лесная мышь 
(бу’аетиз (А.) 
5угайси5) 





*Белозубка-малютка 
(Бипсиу егибси) 
присутствует прак- 
тически во всех жи- 
вых помещениях, в 
приусадебных садах 
и огородах Варде- 


Усатая ночница 


(Муой5 туяастиб) 


Серый хомяк 
(Стсеши5 
пиетаю!5) 





Серый ушан (Р/есо1и5 
аи таси5) 


Малоазиатский хо- 
мяк (Мезостсеи5 
Бтапай) 





Бурый ушан (Р/есо!и5 


Персидкская пещан- 





(РиртягеЦи5 
рилягеЦи$) 


нисского района. аигйи5) ка (Мептопеу (Р.) 
регяси5) 
Нетопырь-карлик Обыкновенная 


полевка (М@сгои5 
(М.) агтаП5) 











Кожановидный нето- 
пырь (РрётеЦи5 
5апи) 





Кустарниковая 
полёвка (М@сгои5 
(Р.) таро") 
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Представители отряда хищных (Сагиога) млекопитающих, нап- 
ример, ласка (Миле а пг’аП5), в поисках мышевидных грызунов не 
только посещает строения человека, но иногда зимует и выводит там 
потомство. Помимо деревень и посёлков, были случаи, когда ласку 
наблюдали в черте городов. Например, в г. Ереване зафиксировалось 
их присутствие в ущелье Норкского массива, в Аване и Зейтуне. 

В ходе исследований, кроме серой крысы (Ваз погуее1си$) и 
домовой мыши (Маз Мизси1$), которые ожидаемо распространены в 
домах повсеместно, были зафиксированы также представители рода 
лесных и полевых мышей (Аро4етиз), из которых чаще всех встреча- 
лась лесная мышь (Ародетиз ига!епз15). Из полевок (М1сгой1$) самой 
многочисленной и распростроненной оказалась обыкновенная полевка 
(М1сго$ агуа|$). В погребах и мусоропроводах в основном мы столк- 
нулись с серой крысой, мыши же обычно встречались в домах и кла- 
довках. В садах больше всех были замечены лесные мыши (Аро4детиз 
ига!еп$1$). Сезонными синантропами оказались полевки, которые удач- 
но переселялись ближе к домам во время урожая пшеницы и других 
злаков. 

В чертах города Раздан было зафиксировано присутствие камен- 
ной куницы (Матез Гота) (Рис. 1), которая охотилась за мышью (Рис. 
2). Кадры были получены с помощью фотоловушек. 

Сначала устройством была зафиксирована мышь (Рис. 2), а потом 
куница (Рис. 1). Последняя оказалась вблизи человеческого дома, го- 
няясь за мышью. Такое «сужение» ареала обитания хищных приводит 
к синантропизации некоторых их популяций. Однако их присутствие в 
чертах города вызывает большую тревогу, так как некоторые хищные 


представители могут представлять опасность и для самого человека. 
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Рис. 2. 
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Оценивая и сравнивая данные исследований по экологии, поведе- 
ния, расселения и распределения синантропных млекопитающих в на- 
ших республиках, начиная с первой половины ХХ столетия (Даль, 
1954) в интересующем нас вышеописанным разрезе, мы можем 
сделать вывод, что за последнее время проявляется заметная тенден- 
ция к сужению распределения мелких и средних размеров хищников 


вокруг человеческих поселений. 
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